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4. Abstract 

Sapp rasa Iona due to acouatlc treatment In tha annular exhaust duct of a model fan hara baan 
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Modal distribution* of tha fan noise source (fan-stator Interaction only) were measured using In- 
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acouatlc Impedance of tha slngle-dagrea-of-freedom treatment was measured In the presence of 
gracing flow* The measured values of mode distribution of the fen noise source, the flow velocity 
profile and tha acouatlc Impedance of the treatment In the duct were used as Input to the prediction 
program. Tha pradlctad suppressions, under the assumption of uniform flow In tha duct, compared 
well with the suppressions measured in the duct for all teat conditions* 

The interact Ion modes generated by the rotor-stator interaction spanned e cut-off ratio range 
from nearly 1 to 7. 
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Except when otherwise defined in the text, the following nomenclature is used 
in this report. 




complex node amplitude 
speed of sound 

core thickness of the acoustic treataent 
acoustic energy 
frequency 
duet height 

transfer function of two acoustic signals at locations x and y 

- /T 

unit matrix 

acoustic intensity 

free space wave number 

radial wave number 

circumferential wave number 

axial wave number (propagation constant) 

segment length 

spinning mode order 

mean flow Mach number 

radial mode order 

mniber of segments 

acoustic pressure 

radial eigenfunction 

radial coordinate 

resistance of treataent 

temperature 

axial acoustic velocity 
reactance of treataent 
axial coordinate 

admittance interface reflection matrix 
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IT J admittance interface transmission matrix 
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KSX Index n in tba bounda r y layer velocity profile. Thus, VEX - 7 is 

for tba 177th p ow er boundary lay velocity profile. 

0 angular location in the cylindrical polar coordinate system 

e angular frequency 2tr f/e 

P duct mall acoustic admittance 

if complex dimensionless eigenvalue k^.r^, where r^ is the 

outer well radius 
6 Dirac delta 

C duct wall acoustic impedance 

n Hf/c-H/X • non-dimensional frequency 

m l/k m non-dimensional axial wave number 

X wavelength 

P e characteristic ispedance of medium 
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P density of air in tho duet 
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(J),(E) duct segment designation 
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Sound suppressions duo to acoustically absorptive treataent in the 
annular exhaust duet of e model fan have been theoretically predicted and 
compared with measurements. 

The predictions were based on the aodal analysis of sound propagation 
in a straight annular fan exhaust duet with axially segasnted treataent. 
Measured values of the node distribution of the fan noise source (rotoir-stator 
interaction source only) end of the acoustic impedance of the treated segaent 
in the duct were used as input to the prediction prograa. The predicted 
suppressions obtained with the assumption of uniform flow (no radial or 
circumferential shear in the flow) compared well with the measured in-duct 
suppressions for ell test conditions. 

Heasureaeats of the acoustic nodes were node on a fan test aodel 
installed in an aaecboie chamber at General Electric Company's Corporate 
Kaseareh and Developme n t Center in Schenectady, lew York. The fan exhaust 
duet was annular in cross section with a faub-to-tip ratio of 0.5. It 
consisted of three axial segments. Acoustic mode probes were located in the 
first and the third segments- which were hard walled. Initial tests were 
conducted with e hard well second segaent. Later measurements were made with 
acoustically absorptive inner and outer walls in the second segment. Tbe 
suppression of sound due to the treated segaent was determined from these data 
by computing the difference between tbe aeoustie energy flux at the upstream 
and the downstreaa aeesureaent locations. The quality of the measured aodal 
data was very good. The acoustic field in the duct at the frequency of 
interest was observed to be stationary over the periods of the aodal 
measurements. This ensured accurate measurement of the modal amplitude and 
phase data. Fan s peeds for the tests were chosen so that the nodes generated 
by the fan-stator interaction spanned a cut-off ratio range from just over cne 
to seven. 




Ik* acoustically* absorbing wall itpasti consisted of single-degrea- 
of-freedoa (XDOF) treataent with liaaar eharactariatlea. Tho measurements of 
tho aoraal acooatic impedance of tho treataant war# carried oat ia Gaaaral 
Electric Company's Acooatic Laboratory at Braadale, Ohio. Thaaa included 
aaasureaants under no-flow conditions oa the treataent hardware and under 
grazing flow conditions oa a sample of the treataent. 

The aodal probe data and a users’ guide to the eoaputer programs are 
available as separata publications (sea list of references). 



. 2.o qaassaa 

Acoustically absorptive traatasnfe in aircraft ansina duets is an 
assantial part of tha overall aircraft noise reduction effort. With the 
increased ea^hasis on energy efficiency, it is necessary to maximize the 
effectiveness of the treataent. Duct treatasnts aust be designed to suppress 
specific noise sources in short duets. By aasiaising the treatment 
effectiveness, a reduction in the treataent length required to achieve a given 
overall engine noise level can be realised, this reduction translates into 
reduced fan reverser length required and, consequently, results in lower 
installed engine weight and higher performance consistent with the desired 
goal of energy efficient noise reduction. In order to achieve this goal, 
analytical treataent design aust be asde an integral part of the advanced 
energy efficient noise reduction technology. 

Recognizing the importance of spinning nodes in the design of the 
ecoustie treataent, several studies have been conducted to. investigate the 
propagation of acoustic nodes in axi-syanetric duets. These studies, however, 
are limited to semi- inf inite ducts with unifora wall impedance. Zoruaski 
(Reference 1) developed an analysis to account for the axial impedance 
changes. This analysis also included both the upstreaa and the downstream . 
travelling nodes. Kraft, et. al. , (Reference 2) measured the acoustic nodal 
distribution in a laboratory nodal inlet duct and designed and tested a 
treataent to verify Zoruaski's analysis. 

This report contains the results of a study on the propagation of 
acoustic nodes in an annular exh au s t duet. It involved the devel op m en t of a 
set of coaputer prograas based on the nodal analysis (Reference 3) of sound 
propagation to predict in-duet suppression due to axially sagnented acoustic 
treataent. It also involved an experimental program designed to determine 
whether the analysis capability was sufficient to describe the physics in 
turbofan aft duct suppressors. Although simplified, the experiment was 
designed to contain the first-order physical effects which were thought to 
govern the far field radiation attenuation due to an exhaust duct suppressor. 




if the experiment vu intended to verify tha sound propagation theory in a 
raalistie turbofan axhaust duet, it contained tha following features: 

1. A raalistie fan stage (RASA Lewis nodal fan— rotor 55) eapable of 
producing well-defined spinning nodes. 

2. An axhaust nozzle with a raalistie eontour and contraction ratio. 

3. A straight annular section axhaust duet to nodal tha simplified 
sound propagation analysis. 

4. A uniform, linear treatment design to nininize impedance 
sensitivity to the operating conditions of tha text vehicle. 

Over tha three fan speeds tasted, spinning sodas ware generated with 
cut-off ratios ranging between just over 1 to 7. The acoustic naasurenents 
involved the neasurenent of (a) modal coefficients upstream and dovnstreaa of 
a treated segnent in the exhaust duct, (b) the aft radiated acoustic far 
field, and (c) the acoustic impedance of the treatment used in the aft duct. 

The test vehicle was nounted in an anechoie chamber in order to neasure 
the far field radiated noise under free field conditions. 

The measurenent of the acoustic inpedanee involved the use of the in- 
si tu nethod (References 4, 5, 6) and of the Acoustic Plunker (Reference 7), a 
non-destructive portable transducer developed by the General Electric Company 
to neasure the normal acoustic inpedanee of finished treatment panels. A 
total of 90 measurements on both the inner and the outer treatment sections 
were made with the Plunker to establish that the normal acoustic inpedanee of 
the treatment was uniform over the area of these sections and that both the 
inner and the outer surfaces had equal inpedanee within prescribed 
manufacturing tolerances. 



The theory of sound propagation in a uniform duet with fluid flow and 
its extension to duets with axial impedance segmentation is described in 
Section 4.1. Section 4.2 deals with the theoretical aspects of the 
computation of the nodal coefficients from the node probe data and the 
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theoretically determined mode eigenfunctions and axial wits numbers in tbs 
bardwall segment of tbs duet. Ths problem has been considered for the unifora 
flow and ths radially sheared flow easss. It involves' ths solution of a sot 
of linear simultaneous equations for the complex nodal coefficients. In ths 
ease of unifora flow the eigenfunction of a downstream propagating (m,n) mode 
is identical to that of ths upstreaa propagating mode of the same order. This 
allows decoupling of ths problsa which ashes the solution siapler in that 
saaller size matrix equations have to be solved. In ths ease of radially 
sheared flow, the eigenfunctions of the upstreaa and ths downstream 
propagating modes of the same (a,n) order are not identical which means that a 
much bigger size matrix equation has to be handled. As the sensitivity of 
matrix solutions to small input errors increases with the size of the aatrix, 
redundant measurements were included in the analysis giving a "least squares 
fit" type of solution from the available data. 

I 


The theory of the in situ iapedancs measurement technique is described 
in Section 4.3. This method assumes that the treatment is point reacting, and 
there is no transmission of sound through the walls of the adjacent cavities. 
The sound field inside the cavity is assumed to consist of plane waves only. 
The method employs the complex ratio of two acoustic originals — one at the 
surface of the treatment and the other at a known location in the cavity 
(usually at the hard back wall). For this reason it has often been referred 
to as the '2-microphone method.* An error analysis of the technique is 
included in Appendix B. 


The description of the experimental apparatus and a selection of 
representative test data are included in Section 5. The test vehicle was 
mounted in the Aero acoustic Anechoic Facility in the General Electric 
Corporate Research and Development Center, Schenectady, Sew York. The airflow 
and noise source consisted of the 13 bladed RASA Lewis 0.504 a diameter fan 
designated as Rotor 55. From the original eleven outlet guide vane set, eight 
vanes were used at a spacing of 0.5 rotor tip chord length from the rotor to 
generate a rotor stator interaction tone at levels appreciably higher then 
noise produced by other mechanisms. To prevent rotor turbulence noise 
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generation, • turbulence control structure was .incorporated, in the inlet duct 
of the fan. The annular flow path of the fan exhaust duet consisted of two 
hardwall sections where in-duct acoustic probes were located, a treated 
section end a noszle termination which was representative of typical engine 
nozzles. The acoustic measurement in the duet involved circumferential 
traversing of an array of pressure transducers through 360 degrees. Signals 
were recorded at 18 degree intervals and then analyzed. In order to ensure 
stationarity of the acoustic field over the tim e of recording the signals from 
the two arrays (of twelve tranducers each) used in the measurement, the fan 
speed was maintained at a constant value during a test. 

Aerodynamic measurements were made to determine the velocity profile at 
each of the in-duet probe locations. This information was necessary for the 
calculation of the axial wave numbers and radial mode shapes (eigenfunc- 
tions). Instrumentation included: three total pressure and total 

temperature rakes, each with five radial stations. In addition, two pitot 
tubes which could traverse radially and circumferentially and aixteen static 
pressure taps mounted in both inner and outer walls were employed.' All 
aerodynamic probes and rakes were removed during acoustic tests. 

The effects of grazing flow and of sound intensity on the Impedance of 
the treatment were investigated using laboratory samples constructed from the 
same materials used in the construction of the annular duet hardware. Steady 
flow resistance measurements were made first after resnving the flexcore (but 
leaving the bonding agent intact). The sample was then instrumented for in- 
situ impedance measurements in the Crazing Plow Duet using a thick walled 
cylindrical cavity to ensure local reaction aspect of the experiment. 

The theory-experiment check is fully discussed in Section 6, and the 
major conclusions drawn from this work are listed in Section 7. The exhaust 
duet suppression prediction program (Reference 8) based on the modal analysis 
was used to predict the suppression due to the treatment in the exhaust duet 
at three fan speeds. The duct was modelled as a three segment straight 
annular duet with the treated segment placed between the two hardwall 
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itsMoti. Mod* coefficients based on ths assumption of uniform flow through 
tbs duet were first obtained from ths in-duet measurements. The mode 
coefficients of the forward (downstream) traveling modes at the upstream 
measurement plane (source plane) and of the backward (upstream) traveling 
modes at the downstream measurement plane (termination plane) were used to 
specify the source and the termination matrices. - With these parameters and 
the knowledge of the impedance of the treated segment, the duct geometry, the 
flow conditions and the spinning mode order number, the program (Heference 8) 
calculates the eigenvalues, the axial propagation constants (both forward and 
backward) , the uniform section transmission matrices for eaeh segment and the 
reflection and transmission matrices of the segment interfaces. The program 
then sets up the stacked system matrix equation and solves it to obtain the 
forward and backward complex mode coefficients in each segment and the modal 
energy flux at eaeh plane. The net energy flux at each plana and the overall 
sound power level (PWL) suppression are then calculated. Suppressions 
predicted in this manner agreed wall with the suppressions obtained from the 
in-duet measurements . 

The sensitivity of the predicted suppression to the treatment iopedance 
values was also examined by varying both the resistance and the reactance of 
the input to the prediction program. The predicted suppression was found to 
be more sensitive to variations in reactance, particularly at the maximum fan 
speed tasted (1900 Hz) when a variation of 0.1 pc. in reactance resulted in 
as much as 4 dB change in suppression. 

Predictions were also made using the sheared flow mods coefficients. 
These did not compare as well with the measured values as those obtained with 
the uniform flow mode coefficients. This is considered to be due to a 
significant energy mismatch that occurs at the segment interfaces when using 
the sheared flow option in the prediction routines. This problem has not been 
resolved. 

Based on the theory experiment check presented in Section 6, it can be 
concluded that the modal analysis for sound propagation in segmented duets can 


be used to predict the in-duct suppression duo to acoustically absorptive 
treatment in the exhaust duet of a turbofan. The aodal distribution of the 
acoustic source and the reflection characteristics of the duet termination are 
required for the prediction. They sap be obtained frost aaasureaents or from 
analytical methods . 


The acoustic data from the mode probes in terns of the linear amplitude 
and phase measured at defined (z, r, 6) locations in the hardwall segments 
of the aft duet are included in the data report (Befereoee 9). 
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3.0 OBJECTIVES 

The objectives of the work reported in this report ere: 

j 

1. Development of the analytical tools required for rapid and coat 
effective evaluation of acoustic treataent designs for fan exhaust 

ducts using nodal propagation analysis. j 

2. Provide experiantal substantiation of the above analytical tools 
in a realistic turbofan exhaust duet eaploying a model fan stage 
and a realistic exhaust nozzle. 

3. Study the effect of the selective reflection by the nozzle on the 
far field suppression. 



ii 
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In this section of ths report, the theoretical model of the aeoustie 
save propagation in an acoustically lined cylindrical duet of unifora cross 
section carrying radially sheared flow and its extension to axial segmentation 
of the treatment is described. - This theory is the basis of the computer 
programs described in Reference 8. In addition, the theoretical bases of the 
computations used in the determination of the modal coefficients from the mode 
probe data obtained in the hardwall segments of the fan exhaust duct and of 
the in-situ impedance measurement technique are described. 

4.1 WAVE-FIELD THEORY 

The prediction of suppression of sound due to treatment in the exhaust 
duet is based on a theoretical analysis of sound propagation in axially 
segmented annular ducts (Reference 3). Propagation in the duct is considered 
in terms of the treatment and the presence of both hardwall and treated 
segments in the duct are incorporated in the analysis by considering the duct 
to be axially segmented and by allowing propagation of modes in both forward 
and backward directions. Reflection and redistribution of acoustic energy at 
segment interfaces and at the exhaust nozzle are also considered. The theory 
of sound propagation in a uniform duct with flow and its extension to ducts 
with axial impedance segmentation are described in the following paragraphs. 

4.1.1 SOUHD PROPAGATION IN ANNULAR DUCTS WITH FLOW 


The propagation cf sound in a duct carrying uniform flow is governed by 
the convected wave equation 


V 2 ? 




(4.1) 


where V ( is the mean flow velocity in the axial direction, V 2 is the 
Laplaelan operator and p is the acoustic pressure. To solve this equation by 
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the separation of variables! as sum that ths acoustic pressure can be 
represented in cylindrical coordinates (Figure 4.1) as 


p(r,e,z,t) - P r < r >P e(e ) # l**”*“t (4.2) 

This reduces the wave equation to two ordinary differential equations, namely 



~*e p ® 


(4.3) 


and 



r dr 



(4.4) 


with 


**(1 - M 2 ) +• 2MkX + (k 2 -k 2 ) > 0 (4.5) 

k is the wave number in free space and k ,k., 'and X are the wave numbers 

r o 

in the radial, circumferential , and axial directions respectively. H is the 
man flow Haeh number (equal to V^/c) . 

Equation (4.3) has a solution of the fora 

P e («) - C e e 1 "® (4.6) 

where values of k^ are restricted to integers since the coordinate 6 is 
periodic with period 2« and the pressure aust be single valued. Physically 
the integral values of k fl , to be denoted by a, represent eireuaferentlal 
mode orders. If the duct contains radial splitters that are lined, k. may 
be coop lex. 
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equation (4.4) with k ft -a is ths standard Bessel aquation and has 
tha solution 


P r < r) - V k r r> + C m V k r P> 


(4.7) 


where and ara Bassal functions (of ordar ■) of tha first and second 
kind respectively. Tha acoustic admittance boundary condition that aquation 
(4 must satisfy at tha duet walls can be expressed in terms of tha 
ci -.«• Inuity of particle displacement as 



where and 0^ are the specific acoustic admittance of the inner wall at 

radius r^, and the outer wall at radius r^. Substitution of equation 

(4.7) into equations (4.8) and (4.9) yields a set of two simultaneous 

transcendental complex equations Which must be solved for k p and C . For 

a jlven circumferential mode order m, these equations have a sequence of roots 

(1c r„) which represent the eigenvalues for a sequence of radial modes. 

The functions p (r) obtained from equation (4.7) represent the corresponding 

eigenfunctions. Setting k r_ » y, r. /r, » p and eliminating c 

r 2 I* a i 

from the simultaneous equations obtained by the substitution of equation (4.7) 
into equations (4.8) and (4.9), we get a single equation for the eigenvalue of 
the form F (?) - 0 (see Appendix C) . 

The eigenvalues can be obtained by solving this equation using the 
second order lewton-Baphson iteration formula 
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J L+r 


T i" 


2F(T i >T' 

21?*^) I 2 - r^Jr-fy) 


<4.l0) 


where y '- v -he initial (guessed) value of the root and y^ + 1 is the 
iterated wiue of the root. The primes denote differentiation with respect to 
the argument. This formula is repeatedly applied until successive iterations 
give roots (eigenvalues) whose absolute values differ by less than 10~^. 

The accurate and reliable determination of these eigenvalues is a 
critical part of the calculation of sound propagation in duets. The ability 
of the solution procedure to converge to the correct sequence of eigenvalues 
depends critically upon the initial value of each root at the start of the 
iteration process. In order to assure convergence, the iteration is performed 
in several steps as indicated below. 

Hectangular duet hardwall eigenvalues are used as starting values and 
the radius ratio is slowly decreased from unity to the annular duct radius 
ratio in order to obtain annular duct hardwall eigenvalues. The latter are 
used as starting values to obtain the eigenvalues for the case of hard outer 
wall and Inner wall of admittance 3^ This is done by slowly incrementing 
the inner wall adadttance magnitude from zero to |0^| and iterating along 
the line of constant phase of 3^. Using the new eigenvalues as the 
starting values and slowly Incrementing the outer wall admittance magnitude 
from zero to 1P 2 I (along the line of constant phase of P 2 ) the 
eigenvalue for the annular duct with inner wall admittance equal to 3^ and 
outer wall admittance equal to P 2 are determined. 

When the inner and outer wall admittances are equal, the two step 
process described above to obtain softwall eigenvalues from hardwall 
eigenvalues reduces to a single step process. In this ease the inner and 
outer wall admittance magnitudes are incremented simultaneously and equally in 
small steps from 0 to |P tf | where 3 W is the adadttance of both walls. 
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The itat iteration proetdur* is ussd while the flow Bach number is 
increaaed in nail stops from xoro to M in order to obtain the lined annular 
duct eigenvalues in the presence of unifora flow. 


This careful step-by-step iteration process provides a reasonably 
reliable eigenvalue routine. However, the acthod is not totally infallible 
and occasionally cases of missed nodes or nodes found more than once are 
encountered. 

An alternate aethod to determine the eigenvalues is used for caaes for 
which the iteration process is not successful. In this aethod the eigenvalue 
equation is set up as a differential equation of the fora (see Appendix C for 
details of the eigenvalue equation) 


F(a,hr 2 ,H,y,pj^ 


(4.11) 


where x is an independent parameter in the equation. To obtain the 
eigenvalues corresponding to wall admittance P^ (for both inner and outer 
walls), equation (4.11) is integrated, using Kunge-Kutta method, from x • 0. 
The step size (Ax) in the integration process is kept small to obtain good 
accuracy, furthermore, in order to minimize the error accumulated in the 
integration process, a second order lewton-Raphson iteration of the eigenvalue 
is performed at each step. Use of this eigenvalue solution pr**cedure in 
combination with the original iterative procedure improves the reliability of 
the eigenvalue solution especially in identifying a soda that could have been 
missed by the iterative procedure. 


After the soft-wall eigenvalues have been determined, the axial mve 
nuaber JZ can be computed from 


1 - (1 


-Mil 


(4.12) 






In the e«s* of m hard wall duet (B^ • B^ ■ 0), whan tha expression under 
tbs radical becomes negative, X bo r oa ms complex and causes the node under 
consideration to decay exponentially. Under such conditions the node cannot 
transport any acoustic energy. Tha cut-on frequency f* for a soda can thus he 
defined by setting 

k* « k^l-M 2 


and obtaining 


f* 


CY Vl-H 2 
2« 2 


(4.13) 


The cut-off ratio Which is the ratio of the nodal frequency to the cut-on 
frequency of the node, is given by 


l 



(4.14) 


where the frequency psramter . k node is considered cut off for 

!<1 but will be propagating when pi. The definition for the cut-off 
ratio as given in aquation (4.14) is valid for a rectangular duct of height 
r 2 if y is replaced by the- rectangular duct eigenvalue It r^, and for a 
cylindrical duet of radius r 2 with y equal to k f r^. 

The effects of sheared flow on the wave propagation are evaluated by 
assisting that the boundary layers at the duct walls are sufficiently thin so 
that the predominant- part of the acoustic energy flux takes place in the 
uniform flow region. In this ease the solution to the differential equation 
for acoustic pressure in the presence of shear layers is approximated by 
choosing the eigenfunctions to be of the same form as in the uniform flow ease 
(i.e. , equation (4.7), with eigenvalues modified by the effects of the 
boundary layers. To obtain the eigenvalues in the presence of wall boundary 
layers, the equation governing the radial variation of acoustic pressure, i.e.. 


* 2 Pr t i»r , 2k /,2 

^2 r dr 1 - H* dr dr y*r “ f 2 J **r “ 


(4.15) 



(where k ■ JC/k » non-dimensional propagation constant) is rewritten as two 
simultaneous first order differential equations, namely 


and 


da. 

dr 



l-M dr 


)q - CkJ 





q 


(4.16) 


and solved numerically by a combination of Bunge-Kutta integration and 
■evton-Kaptason iteration as outlined below. 

I 

l 

The no-flow softwall annular duct eigenvalue (y ■ k r„) and the 

r 2 

corresponding propagation constant e » tfk, as calculated from equation 

(4.12), are used as the initial estimates. At the inner wall p p is assigned 

a value equal to the uniform flow eigenfunction value and the inner wall 

boundary condition equation (4.8) is used to calculate q * dp /dr. These 

r 

values of p and q provide the initial values to integrate equations (4.16) 
r 

by Bunge-Kutta method across the annulus to obtain p p and q at the outer 

wall. The values of p and q at the outer wall are used to calculate tbe 

r 

Admittance 

P - (4.17) 

ik(l-HO p r*r 
C 2 

If this value differs from the specified admittance p^, then the 
process must be repeated with a new value of y until tbe condition 
P>02 is satisfied. This is done by second order Bewton-Baphson 
iteration. 3y setting the difference B-B^ equal to F (y) a revised 
eigenvalue is defined as 


2r(v)r(Y) 

2(F*(y>) 2 - r(y)F"(y) 


(4.18) 
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where F* Ct> - dF/dy end JT <Y> « d 2 F/dy 2 obtained by .« 
performing Bunge-Kutta integration of equations (4.16) with T ±. t and 
Y ± ic as initial aatiaatoa for tba eigenvalue. The Iteration process 
is repeated until an eigenvalue tbat satisfies tbe outer trull adaittance 
boundary condition is obtained. This procedure can be used with any arbitrary 
boundary layer flow profile by appropriately specifying tbe Kaeb number 
gradients dM/dr. Tbe linear or tbe one-sevantb p ow e r law boundary layer 
profiles are used for aost calculations. 


la order to obtain good accuracy of tbe integration process,' a variable 

step slse is used. Logarithmically spaced steps for tbe power law profile 

| 

provide smaller steps in tbe region of bigb velocity gradient. 

Having determined tbe eigenvalues and eigenfunctions, tbe general 
solution for acoustic pressure is constructed in terms of a series 
representation utilizing tbe eigenfunctions, i.e.. 


p(r.O.z) -5n^ A «n p e <#> x 

II t] 

W • “ J 


(4.19) 


vhtr* ths suaastion is dons ovsr all spinnins aodu ordsrs a and radial aode 

ordars n. As discussed before, tha radial variation can ba raprasantad by the 

eigenfunction fora of aquation (4.7). k and K are tha aisanYalua 

r i an an 

(divided by r_) end tbe eziel propegutlon constant respectively for tbe 

* v* 

(m.s) mode. A is designated tbe mode coefficient and is normally 
detersdned by using tbe orthogonality property of tbe eigenfunctions. 
Unfortunately, tbe eigenfunctions in tbe case of wave propagation in duets 
with flow are not orthogonal in the usual Sturm-Llouville sense. A 
generalized definition of orthogonality developed by Craft and Veils 
(Beferenee 10) can, however, be used to evaluate the mode coefficients. 
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4.1.2 MODAL AIALTSXS Kt 8 1CMEM 1 IP DOCT 

The solution for wave propagation in a duet with aultipla traataent 
elements utilises a transfer Matrix principle (References 1 and 2) to eonnect 
the solution at one end of the duet with the other. The duet is as timed to be 
coaposed of axially uniform sections which adjoin at planes where a 
discontinuity in wall admittance ocurs (Figure 4.1). Modal reflections and 
redistributions at admittance discontinuity planes and the consequent 
existence of forward and backward traveling waves in the uniform sections are 
taken' into account in the solution. . 


Based on the modal expansion of equation (4.19) the acoustic pressure 
at an arbitrary axial position z in section J (see Figure 4.2) in the presence 
of forward and backward traveling waves is given by 


P W > (r.d.z) - 1 X 


p e < 0 )P r <*r U> r> 


-ilT* J) (z.-z) 


♦A^ J> e “ j P e (6)p r (k p “ (J) r) 


m n 


(4.20) 


When the duct and liner geometries are axisymmetric, the spinning modes are 
not coupled. The analysis can then be simplified by restricting attention to 
one spinning mode order (m), i.e. , 


,°><r 


J* 0 ' 1 ’ 

uC^z-z. .) .... 

. n j_1 p p Oc p {J) r) 

a a 


n 


(4.21) 
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r » r 



r = r, 


Figure 4.2. Plane and Section Dealgnation Convention for Two Adjacent Duct Treatment Segments 
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Miter* th* spinning order subscript a be* been dropped froa A, X tad k for 

V 

brevity, in equations (4.20) end (4.21) the lower ease superscripts are used 
to indicate values at duet admittance discontinuity planes while th* upper 
ease superscripts indicate values in unifora duet sections. Th* plus (+) and 
ainus (-) superscripts indicate forward and backward propagation 
respectively. A completely analogous equation bolds in section X and other 
sections of th* duet. 


In any uniform section of the duet, sine* eaeb node propagates 
independently, tbs node coefficients at tbe end planes of th* section are 
related by the axial propagation constant and tbe section length. Thus, for 
the radial node n. 


and 



(4.22) 


By representing the nod* coefficients for tu* radial modes as eleswnts of a 
column matrix, th* above equations can be written in tbe following matrix form 



(4.23a) 


for forward propagation, and 



(4.23b) 
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the acoustic pressure and the axial component of acoustic velocity to reaain 
continuous scross the discontinuity. The pressure continuity condition 

p(r. *j) - p(r, x^) (4.25) 

is witten, using the modal expansion of pressure, in the matri. form 

{p + «>! 1 + { ,-<•»} 

+ t4 J ‘ > 

where the superscript T represents the transpose of a column matrix. The 
axial acoustic velocity is expressed in terms of the acoustic pressure (using 
the momentum equation) as 


(r) 


p eOc-M t . ) 
o n 


P„ (r) 


(4.27) 


so that the velocity continuity equation 

v z (r.x j ) - v l (r.x k ) 


(4.28) 
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can b« written as 


H 1 NW 
♦ {,-»}* h 

■ H r M l A+(k, l 

+ H 1 M i 4 ' (k) l 

where [ 0 ] is the axial adaittanee matrix with eleaents 


(4.29) 


„+(J) 


,+(J) 


(n,n) 


k - M X 


+(J> 


(4.30) 


After soae lengthy natrix algebra equations (4.26) and (4.29) can be 
solved to obtain 


and 


H-| 

T +(k).+(J)l ) 

[ ]* l 

| + »~(k) 

M 

(4.31) 

W ; 

x -(J).-(k)j j A -(k)j 

l + 

M 

(4.32) 


where IT] and [B] represent the transmission and reflection aatrices for the 
admittance discontinuity plane under consideration. The deteraination of [T] 
and (B) aatrices requires several matrix algebraic operations as shown in the 
following equations 


24 






with 


lT +«),+(j)j „ ih * {k, * +(,c, j” 1 chJ u> ’ +<j, j 


jg+(K>,+(J)j m |g j jg+OC) ,+{J) j _ jq+<K) ,+(J) j j p+(J)j 
XX s 


[Gjl - P r <,> ]IG +<K), “ <,> ] 


-1 


and 



(4.33) 


where p p is the eigenfunction represented in equation (4.7). Similar 
equations hoXd for other transmission and refXeetion matrix elements. Since 
the eigenfunctions are expressed in teras of the BesseX functions of the first 
and second kind. caXcuXation of the reflection and transaission aatrix 
elements requires the evaluation of the integrals of the products' of the 
BesseX functions. 

Equations (4.23a) and (4.23b) . involving the unifora section aatrix 
lU] . along with equations (4.31) and (4.32). involving the transaission aatrix 
[T] and the reflection matrix [R] . when written for each uniform section and 
discontinuity in the duct, are sufficient to relate the forward and backward 
wave solutions at tme end of the duct with those at the other. Use of this 
method to calculate the propagation of waves in a laboratory exhaust duct of 
finite length would require specification of the axial admittances at the end 
planes to establish the transmission and reflection utrices. (The axial 
admittance at the duct termination depends on the method of flanging the duct 
and the location of the boundaries in the free field. It is further 
complicated by the nonuniform exhaust flow field. At the source location, the 
axial admittance must be related to the internal impedance of the source.) An 
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alternate approach is to choose the two and planes at arbitrary stations in 
hardwall sections of the duet upstream and downstream of treated sections and 
specify the overall modal participation (i.e. , the rum of the forward and 
backward t ravelin* waves) at these plsnes. The modal participation can be 
obtained from in-duct modal measurements. 


The equation which relates the forward and backward modal vectors in 
the source plana is 

N M - [«; <s) -- (s) ] • M <‘.«> 

i 

and at the termination place is 

N M ■ M ■ W <4 - 35> 


where 


(R 1 » reflection matrix at source plane 
s 

(Rj] a reflection matrix at termination plane 
[I] a unit matrix 


sod {Q s } and CQj} are the generalized source and termination 

vectors. If, as suggested above, the source and termination planes are chosen 

as planes in hardwall sections of the duct then 

" [R 8 1 - -II] - tHj] {4.36) 

and the source and termination vectors become the sum of the forward and 
backward waves which can be obtained from measurement for input to the 
program. Equations (4.34) and (4.35) along with equations (4.23a), (4.23b) 
for each uniform admittance section of the duct and equations (4.31) and 
(4.32) for each axial admittance discontinuity fora a completely determined 
system. This system of equations is written as a stacked system matrix 
equation 
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[SJ {A} - {Q} 


(4.37) 


where (S) is the stacked system matrix and [Q] is the generalized source 
vector. Equation (4. 37). .is solved by a double-back substitution routine to 
obtain 

(A) - ISJ — 1 {q} (4.38) 

The stacked system matrix is shown in Figure 4.3 for a three segment duct. 

Equations (4.34) and (4.35), along with (4.36), cannot be used for the 
source and termination planes when the duet walls at these planes are lined or 
when the presence of boundary layer is to be considered. Under these 
conditions, the eigenvalues of the forward and backward modes are different 
and the forward and backward mode coefficients cannot, therefore, be added 
simply to obtain the complex pressure field. A separation of forward and 
backward traveling modes at the source and termination planes can, however, be 
obtained experimentally. This will define the mode content incident from the 
source, i.e., the forward traveling mode coefficients at Plane 1, 

{Qg}, and the mode content reflected from the nozzle; i.e., the 
backward traveling mode coefficients at Plane 6, {Q^}. The stacked 
system matrix equation shown in Figure 4.3 can than be modified by setting 
. [Bg] » 0 » [BjJ and substituting {Qg} and {Q^} for 

{Qg} and {Qg}, respectively. Thus, by using the measured forward 
traveling mode distribution at Plane 1; i.e., {Qg} and the backward 
traveling mode distribution on Plane 6; i.e., {Q^}, the stacked 
system matrix (Equation 4.37) can be solved for the forward and backward mode 

r* 

coefficients at all planes. 

Once the modal solution vector {A} is obtainad, the acoustic energy 
(Eg) at the end planes of all segments is obtained from the following 
expression of the axial acoustic intensity due to Cantrell and Hart 
(Reference 11) 
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(4.39) 


I, Cr.*> - (1+M 2 ) !• lp**l ♦ ^ PP* + P o e Mv^v* 


and 


(z) - 2* 


/ V r * 8> 


rdr 


(4.40) 


The energies at the first and the last planes in the segmented duet can be 
used to calculate the total attenuation due to the treatment In the duct. 


4.1.3 VOZZLZ RETLECTIOH CCEFTICIEHT 


If the termination plane of the fan duct analytical model is 
represented by a plane In the hard wall section of the duct right at the 
entrance to the nozzle, the stacked matrix (*) of the stacked systen equation 
(4.37) can be completely determined only if the nozzle reflection aatrix 
CR^] is specified. The natrlx relates the backward propagating' wave 
just upstream of the nozzle (and due to reflection from the nozzle) to the 
forward propagating wave at the same location by the relation. 


U-< T) ) - (A +m ) 


( 4 . 41 ) 


>o theory adequate for the purposes of this study is currently available for 
the prediction of the nozzle reflection coefficients for multimodal 
propagation. 1*^1 can, however, be determined experimentally. 


The experimental determination of the nozzle reflection matrix 1*^1 
can be complicated due to the fact that the axi symmetric nozzle contraction 
causes a scattering of acoustic modes among radial modes of the same spinning 
mode order, i.e.. a mode incident onto the nozzle can be reflected as a 
combination of several modes. The difficulty in determining the elements of 


the reflection matrix in ttw presence of cross-node scattering among aultiplo 
nodes lias in tha determination of the off-diagonal tans. Separating .forward 
and backward nodes upstream of the nozzle gives a system of n equations (for n 

nodes) in tha form of aquation 4.41 in which the reflection coefficient matrix 

2 

ISj] is unknown; i.e. , there are n equations in n unknowns. One 
hypothetical method of determining the nozzle reflection matrix would be to 
generate a radial node in isolation which, upon incidence to the nozzle, is 
scattered as a distribution of radial nodes. The reflection coefficients for 
the incident node, i.e., S ,, S B . . . . ate., can then be 
determined from the forward-backward nodal separation of the complex pressure 

I 

profile measured just upstream of the nozzle. Bach node, in turn, would have 
to be generated individually. It is, however, impossible to generate an 
isolated radial node in a fan-duct system. 

An alternate approach is to use the measured complex pressure profile 
at a station just upstream of the nozzle for several different distributions 
of the incident nodes and separate the forward and backward node 
coefficients. The measurements made for each distribution of incident nodes 
will produce an independent equation, provided that these incident 
distributions are not linearly related. If a total of n nodes are 
participating in the propagation and scattering phenomena, n different 
distributions of incident modes must be generated. This can be achieved by 
introducing axisymmetric phasing afreets in the duet system, i.e., varying 
treatment length, using liners of unequal admittance on opposite duet walls, 
using node 'Scattering obstacles, Helmholtz cavities, etc., to alter the source 
radial node distribution. Tor each distribution of incident 'modes, n 
equations represented by a matrix equation like equation (4.41) is obtained, 
so that for n different distributions of incident nodes n x n equations will 
be obtained which can be solved for the n x n elements of the reflection 
matrix [B^l. To solve these equations, it is required that the incident 
nodal distributions are not equal or linearly related, to avoid numerical 
problems caused by determinants with two proportional rows or columns. 
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Um the nozzle reflection matrix is determined as nsiMttd above. 
Equation (4.41) can ba written as 

|e -(I)..(I> iU .( t ), . 

Equation 4.42 can than ba used to replace Equation 4.3S so that the stacked 

system matrix equation in Flyura 4.3 will have {QjJ-O. Then, with 

(8 ] set to zero, the matrix equation can be solved upon specification of 

the forward traveling source vector {Q ( }. As discussed before, the 

source vector is specified by the awasured separation of forward and backward 

modes just downstream of the fan plane. 

It should again be mentioned that in the case of an annular exhaust 
duct-noxzle system where spinning modes propagate, scattering of modes by the 
nozzle can be assumed to occur among radial modes only if the nozzle retains 
the axisymmetry of t b .e duct. 

4.2 KATHEKATICAL DETAILS OF MODAL DECOMPOSITION 

Starting with the equation that describes the pressure field in a duct 
with flow. 


H N 

P(r. 6. z. t) - ^ £ [A^p r + (r) e 

n-0 


ik*z 


an r 


(4.43) 


consider measurement of the acoustic pressure amplitude and phase (relative to 
a reference signal • * ~ ). The measurement is suds at a numbs r of 

circumferential positions, 0^ » (jf/K), j ■ 0, 1, 2 . . (2 M - 1), 

whore M » H ♦ 1. Applying the discrete Fourier transform [Reference 4] to the 


Measurements , we obtain the order m coefficient of the Measured 
circumferential pressure distribution (for fixed values of r and s) 


|(r ’ * jg 23 p(r, e t , *) a 


i()ahr/H) 


(4.44) 


B * Of i If e e • f 1 H 

• l 

f 

Applying the discrete Fourier transform to equation (4.43) and comparing with 
aquation (4.44) , wa obtain 


4, * i(k + *) 

< A « p r (r> * 
an 


i(-k" r) 


+ A an p r (r) e 


1 - P_(r,r) 


(4.45) 


The radial and axial mode separation will be considered next. When the radial 

mode shapes, p + (r) and p “ (r) are equal, then the radial mode separation can 
an mn 

be obtained by measuring p_(r, z) at (H + 1) values of r, for a fixed z. The 
results are then used to fora the following matrix equation 


where 


H-i 


p r (r l> p r < r l> p r 

bo ml aN 

P r (r 2 ) 

me 


p r <W 

mo 








. Ul* z - tt " « 

i* . “ ♦ »- . " 

so no 


4+ lk ^ Z A 4- '**■* 

A nl* * A ml e 


(4.46c) 


A+ »i* 

A- • 


+ A tf* 


v r i* z) 

W X) 


®» <r n+l* X) 


(4.46d) 


This matrix equation is then solved for the unknown vector, x. vote that the 
value of x is a function of the axial coordinate, z. 

To separate the order (a, n) forward sad backward propagating nodes, 
the circumferential and radial node separation described above is performed at 
two axial locations, z^ and z^. Denoting the components of x 
corresponding to the (n, n) node by X^Cz^) sad 1 (».) . the following 

matrix can be formed for the desired unknowns A* and a” . 

■n an 


£ X - fe 


(4.47a) 


33 



and 



(4.48a) 


(4.48b) 


where 

A* • * 2 - z 1 (4.48c) 

If we examine the uncertainties in the measurement of k* and A - , v« 

. EM 8n 

in that it is sensitive to the factor [1 - cos(1c* ♦ k~ )Az] in the denoai- 

nQ WMX 

nator. In order to avoid this problea, two things can be done. First, careful 
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•pacing can be chosen such thst Kk*^ + k^)Azl doss not equal an integral 
aultiple of 2«. Secondly, redundant measurements can be introduced, which 
vill significantly improve the nodal decosposition accuracy. The procedure for 
using redundant data to effect a least squares matrix solution will now be ex- 
plained. Consider the general matrix equation 

£ i • 2 (4.49) 

Here £ is a square matrix of order q whose elements nay be complex. £ is 
the vector containing q complex unknowns, and a i» the (complex) input data 
vector of order q. for matrices £, Whose determinant is non zero, there is a 

l 

unique solution for £. However, when 0 is determined from experimental 
measurements, small amasurement errors nay lead to significant inaccuracies 
in £. 


To increase the accuracy in the calculation of £, additional 
independent measurements may be added to the system. Assume that s additional 
independent measurements are added to the system. The size of the matrix £ is 
then- (q +• s) by q, and a is of. order (q + s). While £ remains of order q. 

The system is now overdetermined, and no exact solution for £ is possible. 

The approach used is then to define the remainder x* where 


l-EI-2 


(4.50) 


He now search for the value of £ which minimizes the magnitude of x> 
Voting that x is complex, we define the ma gni t u de squared of x by T> 

r* 

ub*r« 


T - T i T* i - 0^ D* lk 5i (* k 


' °i °ik**k - °ij h °*i + Vi 


(4.51) 


where the * denotes the complex conjugate. He now minimize x by setting 
and 3^/3 equal to zero. The resulting equations can be 
written in matrix form as 
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B** ft I ■ E** S (4.52) 


Where T indicates the transpose.' equation (4.52) is the generalisation of the 
least-squares fit for aatrix equations, and ean be applied to either eqaations 
(4.46a) or (4.47a). 

For duets with sheared flow, ^r (r) and ^r^(r) are not 
identical. Under these conditions, the radial and axial aode separation 
cannot be decoupled in the Banner discussed above, but Bust be carried out 
simultaneously. With only two axial locations, this would lead to in verting a 
aatrix of the order 2 (V +• 1) as opposed to order (> + 1) when decoupling was 
possible (1 - q ♦ s). As the sensitivity of aatrix solutions to satall input 
errors increases with an increase of aatrix sise, redundant aaasureaants 
should be used in this ease also. The resulting aatrix equation is 

S* T £ 2 - S* 1 £ (4.53a) 
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(4.53b) 





As before, the * denotes the coop lex conjugate end T iaplies the transpose of 
the aatrix. 


4.3 THE TS-SXTU IHPEPAHCE MEASUREMEHT TCCHBIQUg 


The In-Situ method (References 4, 3, 6) has been used for the 

measurement of the lonnal Acoustic Impedance of locally reacting treatment 

panels. Dean 4 and Kool and Sarin 5 used this technique on single degree of 

6 r 

f random (SDOP) panels while Zandbergen ausured the impedance of 2 degrees 

12 

of freedom (2D0F) s espies using this method. Zandbergen, et. sl. t also 
neuurtd the impedance of the inlet acoustic treatment of a Fokker F28 
aircraft power plant during flight. 

The theory of this method is illustrated here for an SWF liner. The 
sketch below shows a single cavity with an acoustically hard backwall at x • o 
and a porous face sheet at x ■ d. The walls of the cavity are assumed to be 
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rigid so that no transaission of sound between cavities is allowed. Th« width 
of tha cavity is sufficiently saall to allow only plane wave propagation 
between the baekwall and the face sheet. 


A 



The acoustic pressure at any position x in the cavity and frequency f 
is given by 

p(x t f) » P a i<wt + Cos (Tot) (4.54) 

o 

where k is the acoustic wave number w/c, e is the speed of sound in the air 
filling the cavity and the space outside it, u a 2vf and i a /-l. 

The acoustic particle velocity u in the x direction is related to the 
acoustic pressure by the linearized momentum equation: 









K«t + ♦_> 

Where p_(f) * p • , tb« acoustic pressure at the back wall. 

o O 

■It is assuaed that the acoustic particle velocity is continuous across 
the thin porous face sheet. That is 


u*(f) " u I (f) * _l **]!?** P B < f) (4.57) 

Thus the normal iapedsnee at the face sheet is given by 

p.(f) , p A <f) 

(At) - *~T — - -i eosec (kdj ~77T (4.58) 

A -uj(f)pc P B (f) 

where p A < f) is the acoustic signal at the surface of the face sheet just 
outside the cavity. 

As illustrated by Equation (4.58), the In-Situ aethod requires the 
aeasureaent of the complex ratio of two acoustic signals. For this reason, it 
is often referred to as the Two Microphone Method. 

In the presence of grazing flow, the face sheet transducer signal 
p^(f) aay contain flow noise which nay be regarded as randoa in nature and 
uncorrelated with the acoustic signal. In order to remove this contamination 
of the face sheet signal. Equation (4.58) is modified to 


C A (f) • -1 eosec(kd) . H AB (f) 


(4.59) 


where 


H AB <f) “ p A (f,p B** f) / P B {fJ P B * (f> » tln * * v «r*S«<l complex 
transfer function of the two signals at locations A and B. 
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5.0 nPEMMEKTAL PKOGRAM ACT RESULTS 

This section describes the critical experiment conducted to obtain the 
data necessary to test the validity of the theory described* in Section 4.1. 
This experiment consisted of measurement of the necessary information required 
to calculate the suppression achieved by an acoustically treated section in 
the exhaust duet of a fan vehicle. 

The vehicle was designed for a strong blade-vane interaction tone, 
—ting use of Kotor 55 with only 8 outlet guide vanes ^ closely spaced against 
the 15 bladed rotor. An inlet turbulence control structure was incorporated 
into the anecholc-chaaber test facility so as to reduce the production of 
other acoustic modes. 

The treatment section was designed to be as linear as possible in order 
to establish the value of the wall admittance with minima uncertainty. 
Impedance measurements were made on the actual treatment panels used in the 
test vehicle. Laboratory measurements were made to determine quantitatively 
the effects of flow velocity and sound pressure level upon the acoustic 
impedance. 

Acoustic ’’Mode*' probes were incorporated upstream and downstream of the 
treated section (of* the fan exhaust duct) to establish the acoustic modes 
propagating in both directions at two axial cross-sectional planes, one 
upstream and the other downstream of the treated section. These probes were 
designed to measure the modes expected for the vane-bladc interaction as well 
as all other cut-on modes at the tone frequency of the rotor blades. 

Aerodynamic awasurements were ma de at the same axial planes (by the 
same probe actuators) to establish the mean flow profiles, both radially and 
circumferentially, and to determine the boundary layer profiles incident upon 
and leaving the treated sections. 

The data from the acoustic probes were used in conjunction with the 
aerodynamic measurements to establish the composition of the amplitudes and 
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relative phase* of tha acoustic sodas for us* in the analytical prediction of 
tb* acoustic suppression of the treated section, as discussed in the section, 
"Theory-Experiment Comparison” . 


5.1.1 TEST FACILITY 

The vehicle was tested in an anechcic-chaaber 10.(68 at wide by 7.62 a 
long by 3.048 a high, measured fro* the tips of the foea wedges. The wedges, 
.71 a deep polyurethane foea, provide less than ±1 dB standing wave ratio 
down to 200 Hz. Far- field noise aeasureaents were made by an array of twelve, 
6.35mm-diaaeter, far- field microphones (BAX 4135) located on a 5.182 a radius 
arc. The microphones were arranged at 10* intervals froa 20* to 110* relative 
to the fan exhaust. Table 5.1 summarizes the positions relative to the 
nozzle. Calibration was by piston phone, BAX Model 4220, prior to each 
far-field aeasureasnt. 

The far field microphone levels were recorded on a Sangaao Sabre IV 
tape recorder at 152 cm per second (60 ips) tape speed. ' Power spectra were 
generated by processing ths microphone signals through an HP 5451C Fourier 
analyzer system. Power spectrum levels at blade passing frequency were 
integrated for overall aeoustic power levels. Fan operation was stabilized to 
steady state before initiating any data recording. 

The fan was driven in the exhaust mode as shew in Figure £.1. Figure 
5.2 shows Xotor 55 as it was installed in the chamber. The <iir is inducted 
through a vertical stack end turned through 90* as illustrated in Figure 5.3; 
an inlet turbulence control device was included. 

The effectiveness of the 90* turn at the base of the vertical inlet 
stack and the TCS structure was determined. Kiel probe data are shown in 
Figure 5.4. Losses across the TCS and the flow coatcuring were found to be 
small, and varied by less than 0.5% at any location except at tha top of the 
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Table 5.1. Far Field Microphone Distances and Angles 


Relative to Fan Relative te BOszle 


Distance 

Meters 

Angle 

SSS£SSS 


Distance 

Meters 

Angle 

5.18 

20 


3.54 

36 

5.18 

30 


3.73 

44 

5.18 

40 


3.97 

57 

5.18 

50 


4.25 

1 

69 

5.18 

60 


4.56 

80 

5.18 

70 


4.87 

90 

5.18 

80 


5.18 

100 

5.18 

90 


5.49 

109 

5.18 

100 


5.77 

118 

5.18 

110 


6.04 

126 


i 
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Figure 5.2. Installed Test Vehicle In the Exhaust Mode 
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Figure 5.3. Schematic Diagram of the Teat Vehicle Showing Rotor 55 as 
Operated in the Exhaust Mode 
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Figure 5.4. Kiel Probe Data at the Inlet to the 
Axial Flow Mach Number and Pressure 
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duet toward* the outer wall. Mach numbers war* calculated by using static 
pressures that corresponded to tbe weighted mrigt of the two closest static 
wall taps. The flow was uniform at the aeasunaent locations, the only point 
of significant variation is at tbe top of the duet towards the outer wall 
where the flow is about 3% above the seen. 

5.1.2 FAB VEHICLE 

The vehicle used in this program was the KASA Lewis Research Center fan 
Kotor 55 modified for this program to a close-spaced rotor-OGV configuration 
as shown in Figure 5.5. Bather than the original eleven outlet guide vane set 
spaced at 1.5 chord length eight vanes were used at a spacing of 0.5 rotor tip 
chord length. The annular exhaust duet consisted of two hardsell sections 
where in-duet acoustic probes were located, a treated section, and a nozzle 
termination. Test variations included a hardsell section in place of the 
treated section. 

Pertinent design characteristics of the vehicle are suamarized in 
Table 5.2. 


Table 5.2. Kotor 55 Design Characteristics 


Fan Diameter 50.8 cm (20 In.) 

Badius Batio 0.46 

■umber of Fan Blades 15 

Inlet Guide Vanes Bone 

■umber of Stators 8 

Design Tip Speed (100%) I/rt) . . . . 213 a/see (700 fps) 

Design Fan Speed (100% I//0) 8021 rpm 

Stage Pressure Batio 1.16 

Height Flow 27.0 kg/see (59.5 lbo/sec) 
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s.1.3 akrooxbaiixc xsmuHEixmoM 

Velocity profile measurements were made both up*trta and downstream of 
the traataent section bp aeens of two traversing pitot tubas (radial and 
circuaf erentiel traverse), three total-pressure (and total-temperature) pitot 
tubes each with five radial locations, and sixteen veil static taps. Boundary 
layer profiles were date rained at both locations, up and downstream, by Mans 
of two radially traversing b o und ary layer probes in conjunction with the 
sixteen wall statie taps. 

411 aerodynamic instrumentation information concerning type and 
location is summarised in Figure 5.6. 

Pressure data were i.'dividually sampled through two 48 channel 
scan! valve syst e ms . Temperature data were sampled though an HP scanner and 
digital voltMter. The data reduction program included corrections for 
compressible flow, taking into account specific heat variation with 
temperature and humidity. 

5.1.4 ACOUSTIC "BODE** PROBES 

In general, to determine modal coefficients the in-duct sound field 
needs to be determined. The easiest quantity to measure in a duet is the 
acoustic pressure, which for rotor/stator interaction can be written as 

H 

P(r, •. a. t) - 2 


£ • u “' 
n«0 




-i e" a 

e ] e 


HmO-wt) 


(5.1) 
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Figure 5.6, Schematic Diagram of the Fan Exhaust Duct Showing Instrumentation Locations 
for Aerodynamic Measurements 
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this expression is • solution of ths convected WVI equation using separation 

of variables. Hers ths complex constants and Aj^ are the nodal 

coefficients of the forward and backward propagating (a, n> nodes, p (r) 

■n 

and K are the radial node shape functions and axial wave nuabers, with the 
superscripts + and - corresponding to forward and backward propagation, 
respectively* The other variables aro r, 0, x, t, and », which refer to 
radial location, angular location, axial location, tine, and blade passing 
angular frequency, respectively. 

Determination of the eonplsx model coefficients (aaplitude and phase) 
requires that the eosplex pressure (aaplitude and phase) be neasured. In the 
case where priaarily coherent noise is present, as was the ease here, the 
signal processing technique of phase averaging can be used [Reference 13]. 

This signal processing technique selects from a given input signal only that 
portion Which is coherent with a specified reference signal, and the random or 
noncoherent portion of the input signal is ignored (average approaches zero). 

A phase lock amplifier, which uses this technique, can be used to measure 
in-duet acoustic pressure aaplitude and phase. 

In the ease of uniform flow in either cylindrical or annular ducts with 

hard walls, the radial mode shapes are Bessel functions of the radial order n, 

the eircuaferential mode order a, and of argument (k r) . Also the forward 

+ _ r 

and backward propagating values p (r) and p (r) are equal. When radially 
sheared flow is present, then the radial mode shape is also a function of the 
velocity profile, and forward and backward values p + (r) and p~(r) are not 
equal. 

When pressure, radial mode shapes, axial wave nuabers, probe position, 
and blade passing fequency are all known, the modal coefficients are 
determined by the discrete Fourier transform sad matrix manipulation. The 
details of this modal decomposition are discussed in Section 4.2. An 
important limitation of this process is that the use of only two axial 
measurement positions can lead to large measurement inaccuracies [References 
14, 15]. As is shown in Section 4.2, a third axial measurement position 
significantly improves the modal decomposition accuracy. 
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Prcb* Inscription 

Two in-duct mode probe arrays were designed and fabricated so that all 
cut-on circuaferential and radial nodes including forward and backward 
propagating nodes, could be deterained. A sketch of the probe design and 
geonetry is given in Figure 5.7. Bach array consisted of four radially 
located probes, each probe included three axial Kulite pressure transducers, 
in order to provide a redundant neasurenent at each radial location. Two 
probes Mounted to a strut are shown in Figure 5.8. Bach probe had three slots 
containing recessed Kulite transducers (Model LQ-080-5). Each transducer is 
connected by five wires (0.13 nn disaster) routed through the probe shell and 
grooves in the support struts. Both the probe and strut were nade as snail as 
possible based on transducer size and on the acoustic considerations discussed 
below. The Kulite aodel used was chosen not only because of its snail size 
(1.29 an on each side) but also because of its linearity, stability, and 
relatively high sensitivity. Bach array was mounted on a eireunferential 
traversing actuator, as can be seen in Figure 5.9. 


Disturbance of the sound field due to the present... of the strut and 
probe was of primary concern. Probe scattering can alter the transducer 
readings significantly when an acoustic wave is traveling in a direction 
perpendicular to the probe axis. This scattering problem has been previously 
considered [Heference 16). Also, scattering effects from the strut arise both 
from the blocking of the acoustic motion by the body surface and from the mean 
flow nonuniformity due to the strut thickness. This problem was treated using 
the Taylor transformation (Reference 17], which reduces the low Bach number 
problem to an equivalent stationary medium problem. Based on an analysis 
using that approach, the probe diaawter was selected to be 4.76 am and the 
strut thickness 6.35 an. With these values the error in phase measurement was 
expected to be less than 4.5* and in magnitude, less than 1.01. 


Vortex shedding off the probe support strut was considered, and a 
thickness of 6.35 mn was selected to assure that the frequency was at least 
twice the blade-passing-frequency. 
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Figure 5.7. Scheaatlc Diagram Showing the Radial Locatlona of the Mode Probea and the 
Axial Separations Between the Acouatlc Pressure Transducers on Each Axial 
Probe 
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Figure 5.8 


Photograph Showing Two Axial Probes Mounted on a Radial Strut 


Figure 5.9. Photograph of a Complete Set of "Mode Probes" Showing the Circumferential 
Traversing Actuator and Axial Probes at Four Radial Locations Mounted on 
Two Diametrically Opposed Radial Struts 
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The probe and supporting strut wore designed to withstand the stress et 
velocities in exeess of 0.4 Msch nuaber and to withstand regular handling 
during calibration of the probes. In order to neet these requirements the 
machine screws and the probe shell were Bade of tempered 403 SS with a yield 
strength of 150,000 psi. The natural frequencies of the probe and strut were 
deterained to be well out of range of blade passing frequency. 

A schematic of the Beasureaeat system is given in Figure 5.10. Date 
acquisition was controlled by an HP 1000 eoaputer that saapled pressure and 
teeperature data, and stored in-duct acoustic phase and aagnitude 
Baasuraaents. Two aieroprocessors assisted in this task under direction of 
the HP 1000. One microprocessor controlled the actuator and directed a 
multi-channel scanning amplifier (GS 1566) to sample each mode probe 
transducer signal through a phase-lock amplifier (Brookdeal Model ortholoc-SC 
9505). The phase-lock amplifier passes only that portion of the signal which 
is coherent with the blade passing frequency. The blade passing frequency was 
measured by optically sensing the passage of each blade tip to trigger a 
square wave generator. Both the phase and the magnitude of the coherent 
portion of the mode probe signal were determined in this- manner. The second 
microprocessor controlled the operation of the phase-lock amplifier and 
transferred phase quadrant information to the HP 1000. The HP 1000 also 
applied calibration corrections. The acoustic signals were recorded on a 
Sangamo IV tape recorder. 

The mode probes were rotated in 18* increments around the entire 
circumference of the duct, so that data were measured at 20 circumferential 
positions for each test point. After a test point was completely surveyed, 
the data files of phase and magnitude measurements were supplied to the 
computer pro g ram which then calculated the individual modal coefficients, 
using the method for rmdal decomposition described in Section 4. 

The probes went calibrated for both magnitude and phase. The 
calibration included all wiring and support equipment up to the phase lock 
amplifier. The probes were calibrated by mounting each probe in an impedance 
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Figure 5.10. A Schematic Diagram cf the Data Acquisition end Analysis System Used in the 
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tub* in th* tat axial plan* as a previously calibrated 6.35 aa (0.25 in.) 
Bruel and Kjaar modal 4135 aicrophona and than using a speaker at the other 
end of torn tub* to generate a plane wav* at blade passing frequency. Because 
only a plan* wav* was generatedi th* calibration was valid for both level and 
phase. Uith th* speaker set at a fixed level, th* aegnltude of each 
transducer as read with the phase lock aaplifler was amtehed by use of 
variable gain amplifiers for each transducer, lay difference in phase between 
each transducer and th* microphone was recorded and included as a correction 
in the data reduction pro gram . 

The rotary actuator was used to locate the modal probes 
circ umfer entially. The actuator consists of a stationary section and a 
rotating section. Th* stationary section, mounted on th* adjacent ducting, 
provides th* support for th* motor drive, optical locating devices, and 
transducer wire guides. The rotating section, actuated by a motor driven 
chain, contains the probes. As can be seen in Figure 5.9, th* transducer 
wires feed over a rotating disk and wrap around the rotating section of the 
actuator as it moves. A microprocessor for eaeh actuator controls th* motor 
to drive th* rotating section to a predetermined location. The location is 
determined by a rotary optical encoder (Teledyn* Gurley Model 8625) . Although 
th* mi cr opro c essor controls the actuator, a visual check and manual adjustment 
of the angular location of th* Bode Probe * fas carried out. 

5.2 TTOTCU APtO- lCOqgTIC IESVLT3 

This section presents th* results of the vehicle testing, including: 
aerodynamic parameters, measured modal as^litudes, and farfield acoustic 
data. Besults from a es sur saw n ts of th* treatment impedan c e, both statically 
and with gracing sound and flow, are presented in th* following section (5.3). 

5.2.1 SSLSCTZOB OF FAS SPEED POORS 

The fan was operated at three speed points in order to vary th* cut-off 
ratio of the blade-vane- interaction tones. It was varied from a value 


nur unity (near cut-off) to a relatively high maximum value, over seven, so 
as to have a range of values between these two s xtr s — that would provide 
experimental data for a severe test of the validity of the analytic method 
presented in Section 4. The interaction-tone modes were established by the 8 
outlet guide vanes ehosen for the IS bladed rotor: 

m - 15 ± 8k; k - 0, 1, 2.... 

The cut-off ratios (defined as the frequency divided by the hardsell cut-off 
frequency of the mode) for the propagating interaction sedes at the speed 
settings selected are summarized in Table 5.3. 


Table 5.3. Acoustic Modes Generated by Botor 55 With an Eight Vane 
Stator With Associated Cut-Off Frequencies and Batios 


Mode 

Cut- 

Off 

Freq. 

(Hz) 

Hardwall Cut-Off Batios at: 

1900 Hz 
(7600 rpn) 

1500 Hz 
(6000 rpm) 

1000 Hz 
(4000 rpa) 

-1. o 

270 

7.04 

5.56 

3.70 

-1. 1 

1305 

1.46 

1.15 

— 

7. 0 

1695 

1.12 

— 

— 


5.2.2 FAB PEBFOBMABCE 

The operating map for Botor 55 is shown in Figure 5.11a. The data in 
thie plot shows the repeatability of the performance when three different vane 
designs were used under the QCSSE Program (Beferenee 18]. In the current 
p r o g r a m, the aerodynamic data incorporated for use in determining the flow 
field into and from the test section was used to obtain an approximate 
confirmation that the flow from the fan stage was representative of a typical 
fan vehicle. The test section aerodynamic data are presented and discussed in 
the next section (5.2.3). Data points determined from the upstream and the 
downstream locations are shown in Figure 5.11b, for a fan speed of 7600 rpa 




Figure 5«lla« Operating Map for the Rotor 55 







Weight Flow and Pressure Level 
Determined at Upstream Location 
Weight Flow and Pressure Level 
Determined at Downstream Location 








(95%). As noted on the figure, the two locations gsve different values; this 
was because of the radial velocity profiles used (see Figure 5.12) which aade 
the result dependent upon the particular circumferential location of the 
instrumentation. These data were concluded to be sufficient for the purposes 
needed, that is: the modified fan assembly was sufficiently dose to design 

values of pressure-ratio and air-flow-rate to be a representative turbofan 
noise and airflow source. 

As an interesting sidelight, these same type of data in the original 
fan buildup assembly, revealed serious deficiencies that were traced to 
inadvertent misalignment of the outlet-guide-vane angles. This condition was 
detected and corrected, thereby resulting in the successful measurement of the 
aeoustie modal amplitudes and phases as described in the theory-experiment 
comparison (Section 6). 

5.2.3 TEST SSCTIOS AEBODTXAMIC DATA 

Typical radial velocity profiles at the upstream and downstream 
locations are shown in Figure 5.12; these profiles were determined by the 
traversing boundary-layer rakes which were restricted to a single 
circumferential position. There was sufficient difference in the mean values 
of velocities to cause concern about the validity of the data. 

Therefore, circumferential traverses were then also made with the 
results shown in Figures 5.13 and 5.14, at radial insertion depths of .067 
and .102 m, respectively. These data reveal periodic velocity defects around 
the circumference. At the upstream station the peaks and valleys show a 
circumferential offset between the two radii. The defects also shift in the 
circumferential direction as they move downstream. 

Swirling wakes off the outlet guide vanes were the cause. The 
mixing-out of these defects resulted in the more uniform circumferential 
profile at the downstream location. The particular circumferential location 
chosen for the radial profile measurement was responsible for the discrepancy 
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Figure 3.14. Circumferential Variations of the Mean Axial Flow Mach Number at 101.6 mm 
(4.0 inch) From the Outer Wnll; 7600 rpm 












initially of coneara in Figure 5.12. Finer reaolution of the two 
circumferential traverses is shown in Figure 5.15 Which reveals that there 
were no smaller scale defects. 

The circumferential velocity non-uniformity is of concern because the 
theory of Section 4 only accounts for radial velocity gradients. Averaged 
Mach numbers were used for calculation of axial wave numbers and radial mode 
shapes. 


Most aerodynamic mear^raaents were made at 7600 rpm, although enough at 
6700 rpm and 2800 rpm were sa.si. to assure that the fan performance matched 
previous tests (Figure 5.11b). Figures 5.16 and 5.17 show the similarity of 
the profiles both upstream and downstream at 7600 rpm and 6700 rpm. 

During acoustic runs, 'all probes were removed from the airstreaa. Wall 
static taps and recessed boundary-layer probes were monitored to assure that 
the flow field was the same during the acoustic measurements as with the more 
thorough aerodynamic surveys. Figure 5.18 shows a typical comparison between 
the results from the boundary layer probe and earlier results with a 
traversing pitot tube. 

5.2.4 ACOUSTIC DATA 

This section presents the results of the acoustic modal measurements in 
the test section, at both upstream and downstream locations; measurements were 
made at each of the three fan speeds for the test section with the treatment 
in place and with the treatment replaced by a hardwall (zero admittance) . 

Model Measurement lesults 

For the hardwall case, the results from the two separate mode probes 
(upstream and downstream location) provided a quantitative measure of the 
effect of the non-uniformity of the circumferential flow profile (e.g., per 
Figures 5.13 and 5.14) and of the nozzle reflections upon the resolution of 
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the modal coefficients. A typical result of this cheek is summarized in Table 
5.4 in teens of the relative forward propagating energy (dB). The table shows 
two cases, one considering only cut-on nodes in the nodal coefficient 
determination and the other considering, in addition, cut-off nodes up through 
n » 9 circumferential nodes (maximum obtainable fron 20 circumferential 
measurement locations) and a ■ 3 radial modes (maximum obtainable from four 
radial measurement locations). 

Table 5.4. Comparison of Computed Propagating Energy Flux, 

Upstream and Downstream, for Hardwall Test Section 

Blade Relative Energy Flux* (dB) 


Fan 

Speed 

(REM) 

Passing 

Freq. 

(Hzl 

node 

Index 

Cut-off 

Ratio 

Allowing 

Cut-on 

Only 

Allowing 
Cut-on and 
Cut-off 

7600 

1900 

-1.0 

7.04 

1.2 

2.1 



-1.1 

1.46 

0.8 

0.1 



7.0 

1.12 

0.1 

0.5 

6000 

1500 

-1.0 

5.56 

0.4 

0.2 



-1.1 

1.15 

0.7 

0.3 

4000 

1000 

-1.0 

3.70 

2.3 

3.4 


^Relative energy flux is the difference between energy flux measured 
upstream and that measured downstream. 

The agreement, in general, is very good at the two higher fan speeds, 
and within about 3 dB at the lowest speed; the agreement is best for the modes 
nearest cut-off. (Further discussion of the effects of flow on acoustic 
pressure profiles ratio model components is given in the theory/experiaent 
comparison. Section 6). 

The results of the modal decomposition (using :aximot values of a and 
n of 9 and 3, respectively, and eigenvalues for uniform flow with thin 
boundary layer) are presented in Figures 5.19 through 5.24. The cases are 
listed in Table 5.5. 
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Table 5.5. Index of Figures Shoving Modal Amplitudes and Phases 




Figure 

Huaber 


Test Section 


Modal 

Definition 


5.19a 1900 Bard vail Aaplltude 

5.19b 1900 Bardvall Phase 

5.20a 1900 Treated Aaplltude 

5.20b 1900 Treated Phase 


5.21a 1500 Bardvall Aaplltude 

5.21b 1500 Bardvall Phase 

5.22a 1500 Treated Aaplltude 

5.22b 1500 Treated Phase 


5.23a 

1000 

5.23b 

1000 

5.24a 

1000 

5.24b 

1000 


Bardvall Aaplltude 

Bardvall Phase 

Treated Aaplltude 

Treated Phase 
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Figure 5. 19a. Modal Coefficient Amplitudes Measured at the Upstream and the Downstream 
Locations In the Hardwall Configuration; 7600 rpm. 1900 Hz 
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Flgura 3.20a. Modal Coefficient Aaplltudaa at tha Upstream and Down straam Locations 
In tha Traated Configuration; 7600 rpm, 1900 Hz 
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Figure 3.20b. Phase Value* Corresponding to the Aaplltudea of Model Coefficients In 
Figure 5.20*; 7600 rpra, 1900 Ht 





Magnitude (linear) 



‘Vi 




|Sy»; 
I'fcE •>; 

4 '" : . 

fe 


,$g$ 
lit « , 
Js$ 4 H 

KrC-sfa 

I 

li iV -«V ^ ' 
■*j*V 

■ * 1 

■ . * f. 1 •■-»■ ...■ 

■ l Al '. ■-' 

| pAy \r W 
|yW|®s | 

: * , 1 
V.iKr^f i 

1 *%' -v^i i^h** 

I M,*’’ 



Phase (deg) 


Magnitude (linear) 



O O O O OrOO OOOO OOOO 0»*00 O O O O 

9* N* rt rf r-V*oV* Ki (i N 0) of K irf rf rVdn. rt rf K rf 

I I I I I I I I I I I I 


upstream downstream 

Mode Numbers m,n 

Figure 5,22ft. Modal Coefficient Amplitudes Measured at the Upstream and the Downstream 
Locations In the Treated Configuration; 6000 rpm t 1500 Hz 
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Figure 3.23*. Modal Coefficient Amplitudes Measured In the Upstream and the Downstream 
Locations In the llirdwall Conf Ijuratlon; 4000 rpm, 1000 lit 
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Piguro 5.24a. Modal Coefficient Amplitudes Measured In the Upstream and the Downstream 
Locations In the Treated Configuration; 4000 rpm, 1000 H* 
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Figure 5. 24b. Phase Values Corresponding to the Modal Coefficients in Figure 5.24a; 
4000 rpm, 1000 Ha 




The figure* show the modes expected from rotor/stator interaction, 
although other nodes are present also. For the 1500 Hz case only the (-l v 0) 
and (-1, 1) should be present. Figure 5.21a shorn these nodes were indeed 
measured but there was also a strong (-5, 0) node. This would be possible if 
the inlet flow were not uniform, and nay in fact be the result of interaction 
between the one high velocity zone near the top of the duct (Figure 5.4) and 
the fan. The fact that this node is so strong is due to the fact it is just 
cut-on. Also because it is just cut-on it exhibits large reflections, which 
nay explain its drop in level from upstream to downstream. In other cases the 
expected rotor/stator modes were stronger. For example, the forward 
propagating (7, 0) node was about 7 dB stronger than any other mode in the 
19C0 Hz hardwall case. 

Figure 5.22a (1500 Hz) is a typical case with treatment. As expected, 
the upstream forward propagating coefficients are about the sane as in the 
hardwall case, but the downstream coefficients are lower, due to the 
treatment. Figures 5.20 and 5.24 show similar results for the 1900 Hz and 
1000 Hz cases. As Is evident, the treatment proved to be optimally designed 
for 1900 Hz. As expected, in the 1900 Hz case the (7, 0) mode which is just 
cut-on is the dominant mode, although in the 1000 Hz case the (-1, 1) mode 
which also is just cut-on is weak and the (1, 0) mode dominates. 

Treatment was expected to have no effect on the upstream forv»**d 
propagating modal coefficients. Figure 5.25 shows the effects of treatment 
and temperature on these coefficients. Hun A66 was a 7600 rpn hardwall case. 
Run A67 was the same as A66 except for the fact the upstream and downstream 
nodal instrumentation was switched. When run A67 was carrier oun the outside 
temperature was lover and because of this temperature change the (7, 0) node 
is less cut-on for A66 and therefore stronger, while the other modes are in 
relatively good agreement. Hun A71 vas the 7600 rpm treatment run with 
temperatures similar to A67 and it shows good agreement with the hardwall 
cases except the A66 (7, 0) mode as would he expected. It should be noted 
that the A67 downstream coefficients did not make sense and as far as can be 
determined the only possible cause was calibration drift. Run A67 was the 
only run where the probes were not calibrated the day of the run and the 
downstream probes had a different set vi calibration trim potentiometers. 

This could explain why the upstream probes still gave good results. 






In general, the differences between carrying eat tbe nodal 
decomposition using sheared flow values of ^ sad or unifora values 
were saell, although the 1000 Hz case did show a difference of 1.5 dB in the 
(-1, 0) node. The ssae could be said Whether the aodal decomposition was 
solved for cut-on nodes only or for all nodes. 

PaivPiold Measurenent kesults 

Table 5.6 shows the sound power levels for both hardsell and treataent 
cases along with tbe resultant suppression at each naming speed. The 
integration for acoustic power level was done using the nozzle exhaust as the 
acoustic souree and the revised distances and angles fron Table 5.1. 

Table 5.6. Sound Power Levels Proa Par Pield Data 
(Power Spectrua Bandwidth - 4 Hz) 


Blade 

Passing 

Preq. PUL (dB) PUL (dB) Treataent 

- (Hz) acftail TCMfomt Suppression 

1000 120.0 118.8 1.2 

1500 129.2 125.0 4.2 

1900 126.8 116.6 10.2 


Pigures 5.26 through 5.28 show the directivity patterns of the narrow 
band SPLs. The lack of suppression in tbe 30* microphone position could well 
be the result of the directivity shifting. Which would be expected, in 
general, the SPL et blade passing frequency was 5 to 10 dB higher than its 
harmonics for 7600 rpa and 6000 rpa, while for 4000 rpa the difference was at 
least 15 dB. 

5.3 dgggmc-BBAiBari. atm la 

In the theory experiment cheek of Section 6, an accurately determined 
value of the acoustic impedance of the treataent sections was required. 
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Figure 5.26. Plots of the Far Field (Narrow Band) Sound Pressure Level (SPL dB) Against 
Microphone Angular Position Relative to the Nozzle Axis for Hardvali and 
Treated Configurations at 7600 rpm (1900 Hz) 
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Although tho treatment was dosignsd to ho linear, measurements wars carried 
out to ostablish tho offsets of sound intensity variation and grazing flow. 

In addition a comprehensive survey of tho impedance was asdo orsr tho surface 
of both tho inner and tho outer treatment sections to ensure that (a) the 
acoustic impedance of the two surfaces was equal and (b) the treatment 
impedance was uniform over the entire surface. The details of the design of 
the treatment, the measurement techniques and apparatus used and impedance 
data obtained are presented below. 

5.3.1 TWATHOT DSSIGB 

The acoustic treatment was designed to be locally reacting and to have 
equal impedance on the inner and the outer walls of the treated section of the 
aft duet. A single degree of f re edom (SDOF) design was chosen consisting of a 
25.4 am thick Flex Core bonded to a 1.143 am porous faee sheet and the wooden 
back wall of the treatment sections. The Flex Core used is made from 0.063 no 
thick a turn 1m as sheet. The geometry of the calls is illustrated in Figure 
5.29a. 


The faee sheet consisted of a 200 x 400 wire mesh bonded onto a 30% 
porosity aluminum perforate (.83 ns hole size) sheet. The purpose of the wire 
mesh was to make the treatment linear; that is, the effects of flow and sound 
intensity on the impedance of the treatment were expected to be small as a 
consequence of the wire mesh. 

From previous experience of the blockage due to bonding the core to 
face sheets, the resistance of the treatment was expected to be ap pr oxi mately 
1.0 pc (> 41.4 cgs Kayla) . The resistance actually maasured on the 
finished treatment sections and on the laboratory samples was a pp r ox i mately 
0.5 pc. 


The gecamtry of the treated segment of the aft duct is schematically 
shown in Figure 5.29b. 
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Figure 5.29b. Schesatlc Dlagraa Shoving the Acoustic Treatment 

In the Straight Annular Part of the Fan Exhaust Duct 
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s.3.2 xm-sxtu xmfsdahcs «usmvsant 

• HMMfMMt iBomtm 

The affects of grazing flow sad sooad intensity wore studied ia tbs 
Grazing Flow Duct using a saaple of the fees shoot of tho trootaoat. The 
saaple was prepared using tho saae aatorials aad proced u res as used la tho 
aaaufaeture of tho trootaoat sections for tho aft duet, the Flex Coro was 
reaovod while oasuriag that tho blockage of tho face shoot duo to tho bonding 
rosin was preserved oa tho sample. A steady (dc) flow aoasuroaont on tho 
saaplo was aado prior to its preparation for Xn-Situ aeesureasnts. 

Figure 3.30 shows the tost soaplo. The cavity consists of a thick 
brass cylinder with inside diaaoter of 17.78 aa. ' The depth of tho cavity 
could be varied between 23.4 aa and 43.72 aa. Bndeveo pressure transducers 
(Model 8S14-10) were used to aeasure the aeoustie sign els st the face sheet 
and tho back wall. Those transducers are of 1.387 aa diaaoter. The range of 
those is 0-88930 pascals (0-10 pel) with a sensitivity of .004569 ± .0013 
aV/paseal (31.5 ± 9.0 aflT/psi). Tho size and tho sensitivity of thess 
transducers aado the* suitable for the Xn-Situ aeasureaoats. 

Figure 5.31 shows tho data acquisition and analysis hardware used in 
tho aoasuroaonts. Tho signals froa tho transducers were amplified by 
Tektronix aaplif iers (Model AM 502) . Those are DC-coupled differential 
aaplifiors with good coaaon-aode rejection capabilities (rejection ratio of SO 
dB DC to SO kHz) and high gain for low voltage aoasuroaonts. The aaplified 
signals were saapled and analyzed by the Tiae Data Analysis systea based on a 
12 bit* 2 channel A to D converters and a DSC PDF 1135 alnicoaputer. The A to 
D converter of the Tiae Data Analysis systea has aaplifiors aad antialiasing 
filters for each data channel which can be selected through software. This 
ensures that speetral information is uneoataainated by higher frequency signal 
and that the best use of the dyneaie range of the A to D converter is aado. 


Figure 
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Scheaatic Dlagraa of the Apparatus Used 
In the Grazing Flow Duct 


for In-Sltu Iapedance Measurements 
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Pfpurr S.1|, Schematic DlaRr/im of the Data Acquisition and Analysis Hardware Used In the In-Sltu Impedance Measurements 
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In— Situ iapadanea Muarowoti wrt also aada on tho treataent soetions 
prior to their installation in tbo vehicle. Tbo outer traatasnt saetion was 
aada up froa six cireuafarential segaents and tba innar saetion was aada up 
froa tbraa segaents (saa Fignra 5.32). A single In-Si tu asasuraasnt was aada 
in aaeb sagaant. Tba purposs of tbasa aaasuraaants was to astablisb an 
avaragad iapadanea value (without grazing flow) for tba traatasnt to ba tasted 
in tbo aft doet. Tba axpariaantal arrangaaant usad in tbasa asassreaants is 
shown in Figure 5.33. Broadband noisa was usad for tbasa aaasuraaants. The 
aaasuraaants on tba innar saetion wars dona using a saapling rata of 12,800 Hz 
and a bandwidth of 12.5 Hz whila those on tba outar saetion wars carried out 
with a saapling rata of 25,600 Bz and band width of 50 Hz. 

a Imrt d ir e « Data Froa tba Treated Sections of tba Aft Duct - Ho Flow 

Has i stances and reaetanees froa the aaasuraaants on the outar treataent 
saetion are presented in Figure 5.34. The narrowband spectra of the signals 
aaasurad at tbo faea sheet and tba back wall of the cavity at location A are 
shown in Figuro 5.35. 

Tho eohorenca of tho two signals and the phase difference between then 
is shown in Figure 5.36. 

The quality of tba data in Figures 5.35 and 5.36 is typical of these 

aeasuroaents. 

With tbo exception of tho iapadanea aaasurad at cavity B, tbasa 
aaasuraaants look good. They show that tbo resistance varies by a saall 
aaount between tbasa cavities. For axas^la at 1900 Hz, tba resistance 
aaasurad at cavity A is 0.48 pc while that aaasurad at cavity B is .67 pc. 

As tba inner treataent section is aada up of three cireuafarential 
segaants, only three In-Sltu aeasuroaents ware nude on it. The data froa 
these aaasurasMnts Is shown in Figures 5.37. Tba reactances do not vary froa 
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Figure 5.32. Photograph of the Inner Treatment Section of the Acoustically Treated 
Segment of the Aft Duct 
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Figure 5.33. Experimental Arrangement Deed In the In-Situ Iapedence Measureeent 
on the Outer end the Inner Treatment Section# of the Aft Duct 
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Figure 5.34. In-S I tu Impedance Measurements at Six Different Circumferential Locations 
on the Outer Treatment Section of the Aft laict (continued) 
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Figure 5.35. Auto Correlation Spectra of the Acoustic 
Signals in the Tn-Situ Measurement of 
Cavity A of Figure 5.34 
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point to point, bat tho resistance nt cavity c is appreciably higher than nt 
cavities t or b. This observation raised the question as to how 
representative the measurements at a, b and e wore of the iapedanee of the 
panels on which they were located. There was a need to sake several 
additional aeasaresMnts to get a good statistical measure of the averaged 
iapedanee value of the treatment and the standard deviation frew it. To do 
this by the In-Si tu method would be tedious and tiae consuming. Consequently 
it was decided to carry out a comprehensive survey of the treatment impedance 
using the Acoustic Plunker (See S.3.3). It should be noted however that at 
any point the acoustic resistance measured by the In-Si tu method is virtually 
constant between 1000 Ha and 2000 Hz, the range of Interest in this study. 

e Impedance Measurements in the Craxlnz Flow Duct (CTO) 

These measurements were carried out on a laboratory sample of the 
treatment. First the effects of sound intensity were measured without grazing 
flow. Then the effects of grazing flow on the impedance were measured. It 
was assumed that the correlations obtained from these measurements would 
accurately represent the grazing flow and sound intensity effects occurring on 
the treatment in the aft duct during andal measurements. 

e na_ittisyL2L^g3aOat8agiti 

Figure 5.38 thorn tho increase in acoustic resistance relative to the 
value Measured at 125 dfi plotted against the sound pressure level at the three 
frequencies of interest. At 1000 Hz f the effect of increasing SPL on the 
acoustic resistance was such lover than that Measured at 1500 Hz and 1900 Hz. 
Negligible effects of sound intensity changes were observed on the reactance 
of the saaq>le. 

The above data was obtained with discrete frequency sound field in the 

duct. 
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O 1000 Hz 



Figure 5.38. Increase in Acoustic Resistance With Sound Pressure Level (SPL f dB) of the 
Discrete Acoustic Signal; Ho-Flov 
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A cavity depth of 4.572 cm I n ctoe d of 2.54 cm «u ehosea for these 
Measurements because an error analysis (see Appendix B) had shown that the 
Magnitude of the error in the Measured resistance would depend on the 
Magnitude of iapedance at the frequency of interest. By increasing the cavity 
depth* the Magnitude of the reactance of the staple was reduced significantly 
between 1000 Hz and 2000 Hz to yield aore accurate aeasuresMnts of the 
acoustic resistance in this frequency range. It was assuaed that the 
reactance of the cavity did not influence the effects of grazing flow on the 
luap iapedance of the porous face sheet. 

The effects of grazing flow on the resistance of the saaple are shown 
in Figure 5.39a. As expected, the increase in the resistance with increasing 
aean flow Haeh nuaber is roughly the sane at the three frequencies of 
interest. The data correspond to roughly 140 dB (OASPL) at the surface of the 

saaple. 

Figure 5.39b shows the corresponding plots of reactance against the 
man flow Kach nuaber. 

Figure 5.40 shows typical narrowband spectra of the signals sensed by 
the transducers at the face sheet and the back wall of the cavity at Hach 0.4. 

The effect of sound intensity on the iapedance in the presence of 
grazing flow was found to be such less than that Measured under no flow 
conditions. In view of this the sound intensity effects under test conditions 
in the aft duct my be neglected. 

5.3.3 FLUBEER HEASOBEHEBTS OB THE TBEAIHEBT SECTI0E3 OF THE AFT DUCT 

As a considerable aaount of variation In the iapedance of the treataent 
sections was revealed by the In-Si tu Measurements (as illustrated in Section 
5.3.2). A Much larger nuaber of iapedance aeasuremnts on these sections was 
mde using the Acoustic Plunker which is described below. 
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Figure 5.39a. Increase in AcousClc Resistance Vlth the Mean Flow Mach Number in the 
Grazing Flow Duct 



97 





• 1000 Hz 
A 1500 Hz 
.2 r X 1900 Hz 


u 



0 .1 .2 .3 .4 


„ mean FLOW MACH NUMBER, M 
AX m" (X M ' W* S?L-140dB 


Figure 5.39b. Increase In the Acoustic Reactance With the Mean Flow Mach Nuaber In the 
Crazing Flow Duct 
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Figure 5.40. Auto Correlation Spectra of the Face Sheet Signal A 
and Backvall Signal B During an In-Situ Impedance 
Measurement at M * 0.4 


The Acoustic Plunker [Reference 7] is a transducer developed b y the 
General Electric Company for the measurement of the normal acoustic impedance 
of fir* : \\... treatment panels. It is a portable and nondestructive device. It 
has to b«, calibrated to measure the normal impedance of treatment panels of a 
particular design. It has sufficient sensitivity to measure small variations 
of the normal impedance on a given panel at a given OASPL as veil as to sense 
impedance changes with sound intensity. 

The Plunker is schematically illustrated in Figure 5.41a. It consists 
of a tube of circular cross section. The mils of the tube are acoustically 
hard and smooth. A sound source is located at one end of the tube. The other 
end is placed against the surface of the treatment panel such that the tube is 
at right angles to the surface. The rubber flange is designed to minimize 
sound leakage. The diameter of the Plunker tube is such that over the 
frequency range of interest 9 only plane waves can propagate in it. Two 
transducers located at distances and x^ from the treatment surface 
sense the stationary sound field in the Plunker tube. The signals from the 
Plunker transducers are analyzed to yield the time averaged values of the 
complex transfer function: 


H 21 (f) * P 2 (f) P l <f) 7 P l (f) P l (f) (5.2) 

The apparent impedance measured by the Plunker using the above transfer 
function is given by: 


JSln(tot 2 ) - H^fJSinCto^)} 
C p (f) - -i |co*(tac 2 ) - H n (f )Co«Ooc 1 ) } 


(5.3) 


If a sample of the treatment is cut out of the panel at the location where the 
plunker measurement ( p (f) is made and its normal impedance ((f) is 
measured in the Impedance Tube schematically illustrated in Figure 5.41b, then 
the calibration of the Plunker is computed from the above data &3 follows: 


i p (n - 


an 

V f) 


(5.4) 
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The calibration of the Plunkar as defined abova for a givan design of 
treatment is naodad baeausa tha apparent iapedanee C p (f ) measured by tha 
Plunkar is diffarant from tha Zapadanea Tuba data for tha following reasons: 

1. Tha absorber volume behind tha panel surface in tha Plunkar 
measurement is significantly larger than in the equivalent 
Iapedanee Tuba messuresmnt. 

2 . Some sound leakage aay occur due to iaparfact seal be.vjen tha 
Plunkar and tha panel surface due to tha curvature of tha traataent 
surface. Tha affect of this leakage can be calibrated out only if 
tha Plunkar axis is always perfectly noraal to tha surface and that 
tha radius of curvature of tha surface is constant everywhere. 

3. There aay be transmission of sound in the core of tha panel at 
right angles to tha sound field in the tube due to 

(a) non rigid cavity walls 

(b) drainage holes between cavities 

(e) iaparfact bonding of core to tha face sheet and back wall. 

Once the calibration 5p(f ) for a particular panel design is 
obtained in the manner. described above, tha normal impedance of any other 
panel of that design can be measured by simply measuring C p (f) with the 
Plunkar and then multiplying it with 5p(f ) . 

As it was not desirable to cut out a sample from the treatment sections 
for the aft duct, tha Plunkar calibration was carried out using an In-Situ 
measurement of impedance. It was established that this measurement was 
similar to an Impedance lube measurement on a laboratory sample of the 
treatment. The inside diameter of the Plunkar tube used in these measurements 
was 17.78 mm. This size was regarded sufficiently large to give a measurement 
of impedance on a statistically significant area of the panel. It was also 
sufficiently small compared to the radius of curvature of the treatment 
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surfaces which could ho regarded essentially flat over tho area of tba tuba. 

It alto a naur ad mini— leakage of sound hatwaan tha Plunker and fcha faea 
ahaat due to imperfeet saal causad by tha curvature of tha treatment surface. 
Separata calibrations of tha Plunkar were required for tha outer and tha inner 
sections due to significant differences in the curvatures of their surfaces 
which affected tha extent of leakage between tha plunkar and tha treataent 
surface. 


On each of the six segments of tha outer and the three segments of the 
inner treataent sections, at least ten iapedanee measurements were obtained 
with tha Plunkar. During each measurement tha plunkar was simply placed on 
tba location. Apart from its own weight, no additional pressure was applied. 
Care was taken that tba axis of tha plunkar was normal to tha surface of the 
treatment. Because of tha curvature of tha treatment surface the leakage of 
sound was a problem. A very small tilting of tha axis of tha Plunkar from the 
normal causad significant leakage. This problem was particularly severe in 
tba case of tha inner treatment section. This leakage of sound affected the 
accuracy of tha impedance data at lower frequencies (below 1500 Hz). 

Figure 5.42 shows the results of the impedance survey of the outer 
treatment section. It is based on sixty different amasurements. The standard 
deviation of the impedance data is quite small at 1500 Hz and 1900 Hz but at 
1000 Hz the deviation in the data- is much larger due to reasons mentioned 
above. 


Figure 5.43 shows the results of the Plunker measurement on the inner 
section. As expected, the standard deviation of the resistance values 
measured on the inner section is larger than that on the outer section. 
Moreover the resistance values below 1500 Hz are lower than expected due to 
sound leakage. 

Plunker measurements also confirmed that the In-Situ measurement on 
cavity J of the inner treatment section is representative of the impedance of 
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th* entire panel and ie midway between the In-Si tu measurement result* on 
cavities H and J. The Plunker results show lower frequency dependence than 
the Xn-Situ data. 

In the light of the Plunker measurements the following conclusions may 
be drawn. 

1. The averaged impedance of the outer treatment section is 
a pp r oximately equal to that of th* inner section. 

2. The treatment impedance is uniform over th* entire surface within 
the limits: of manufacturing tolerances. 

It was decided that the averaged acoustic impedance values obtained by 
the Plunker at the frequencies of interest should be used in the 
theory-experiment check. These values are presented in Table 5.7. The 
In-Si tu measurements had shown that the resistance of th* treatment was 
essentially constant between 1000 Hz and 2000 Hz. In view of this, th* 
resistance at 1000 Hz was assumed to be the earns as that at 1500 Hz. It 
should also be noted that the Plunker measurements were obtained with an OASPL 
of 140 dB (approximately) at th* treatment surface. As the in-duct sound 
levels during th* tests at Schenectady were significantly lower, appropriate 
corrections to the resistance values would be required.* In addition 
corrections should be made for the effects of grazing flow. The following 
approach is suggested: Using the sound intensity correlation of Figure 5.38 

and the average impedance at 140 dB from Table 5.7, the no-flow resistance at 
130 dB is appro ximatey 0.44 pc. To this should be added th* grazing flow 
effect from Figure 5.39a. The^reaetanee values of Table 5.7 should only be 
corrected for flow effects using Figure 5.39b. 
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, 4.0 tmoBr-ppgmiBrr oomparisor 

The annular duct suppression prediction program baaed on the nodal . 
analysis described in Section. 4.1 was uaed to prediet the suppreaaion due to 
the treatment in the Rotor-55 exhaust duct for ell three test KPHs. The 
theoretical predictions are compared with the naasurenenta in thia section. 
The physical phenomenon of propagation in ■ assented ducts is discussed with 
reference to thia theory-experiment coaperison. 

6.1 JTOT TO PRRDICTIOR PB0CR4M 

The exhaust duct was modeled as a three segment straight annular duet 

I 

with the treated segment placed between two hardwall segments. The lengths of 
the hardwall segments were taken to be the distances of the mode measurement 
planes from the treated segment. The effective length of the treated segment 
was found to be 0.23 meters (L/H *1.8) after accounting for the bole blockage 
at segment ends and treatment panel interfaces. 

The principal input parameters for the program are the duct geometry, 
the mean flow conditions, the acoustic frequency, the inner and outer wall 
impedances (or admittances) for all segments, the spinning mode order, the 
radial mode distribution of the acoustic source and the reflection 
characteristics of the duct termination. The mode coefficients of the forward 
traveling modes measured at the upstream plane (or the "source** plane) and the 
node coefficients of the backward traveling modes measured at the downstream 
plane (or the "termination" plane) were used to specify {Q^} and 
{Qj} respectively. As mentioned in Section 4.1.2 this required that 
(Rgl - 0 * [R^J. For the blads passage frequencies corresponding to the 
rotor RPBs in the test, a maxima of two radial modes (i.e. , n«0 and n*l) were 
expected to be cut-on in the hardwall segment of the duet. Vhile utilizing 
the prediction program, however, a minimum of four radial modes (i.e., n-0, 1, 
2, 3) were considered to participate in the transmission and redistribution of 
acoustic energy in the duct. The acoustic impedance of the treatment in the 
duet was initially assumed to be equal to the average normal impedance 
measured in the absence of flow using the acoustic Plunker (Section 5.3.3). 
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with these input paraneters tbs' proem calculates the eigenvalues, the 
axial propagation constants (for both forward and backward propagating nodes) , 
tha unifom section transaission matrices for aaeh duet segment, and the 
reflection and transaission aatriees of the segaant interfaces. The program 
then sets up tha stacked system matrix equation sad solves it to obtain tha 
forward and backward coaplex node coefficients and aodal energy fluxes at each 
plans. The net energy flux at each plane and tha overall PUL suppression are 
then calculated. 


6.2 






Hode coefficients obtained from the in-duct aeasurenents and based on 
the assumption of uniform assn flow in the duet were first used in the 
predictions. 


Table 6.1 lists the neasured and predicted in-duct suppressions for the 
fan-stator interaction nodes at the three blade passage frequencies 
corresponding to the test conditions. These suppression values represent the 
difference in tha forward energy fluxes between the downstrean and upstream 
neasurenent planes. The effect of reflections froa the nozzle is taken into 
account by specifying the backward traveling node coefficients at the 
downstrean plane. The predicted suppressions can be seen to be in good 
agreenent with the neasured suppressions for all three frequencies and for 
both spinning node orders. The largest difference between the predicted and 
neasured suppressions is 2.7 dB which is cot large compared to the measured 
suppression of 16.8 dB at that condition, lotice that six radial nodes were 
used in the prediction for n_- 7 case (for reasons discussed later in this 
section) while only 4 radial nodes were used for n ■ 0-1. 


The predicted and measured mode coefficient distributions of forward 
propagating acoustic energy at the downstrean plane are conpared in Figures 
6.1 through 6.4 for the four cases listed in Table 6.1. The magnitude as well 
as the phase of the neasured node coefficients can be seen to be in good 


Table 6.1. Comparison of Measured Suppressions with Theoretical Predictions 
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Plyuro 6.1. Measured and Predicted Conplex Mode Amplitudes at the Downstream Plane - 

Forward Propagating Modes - 1000 Hz. m ■ -1. M ■ 0.21. Uniform Plow. Soft Hall 
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PIgurfc 6.2. Measured and Predicted Complex Mode Amplitudes at the Downstream Plane - 

Forward Propagating Nodes ~ 1300 Hi. m • ~1, H * .32. Uniform Plow, Soft Wall 
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Figure 6.3. Measured and Predicted Complex Mode Amplitudes at the Downstream Plane - 

Forward Propagating Modes - 1900 Hz, m m -l, M ■ 0.4, Uniform Flowi Soft Mall 
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Figure 6*4* Measured and Predicted Complex Mode Amplitudes et the Downstream Plane ~ 

Forward Propagating Modes - 1900 He, ra ■ 7 $ M ■ 0,4, Uniform Flow, Soft Wall 
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agreement with the predicted value* for the cut-on aodes. For a » -1 at 
1000 Hz and for ■ • 7 at 1900 Hz the higher order radial aodes (n > 0} are 
theoretically cut-off in the hard wall duct. This causes the predicted 
aaplitudes of these nodes at the dovnstreaa location to be several dB lower 
than the predicted aaplitude of the first radial node. The neasureaents also 
show the aaplitudes of the higher order nodes to be at least 10 dB lower than 
the aaplitude of the cut-on aode. The difference in the aeasured and 
predicted values of the aode coefficients for eut-off nodes will therefor* 
have only a saall effect on the total suppression in the duet. 

Figures $.5 through 6.8 compare the predicted values of the aode 
coefficients for the backward traveling nodes at the upstreaa plane with the 
corresponding aeasured values. The comparison for the backward traveling 
aodes is not as good as that for the forward traveling nodes. For a ■ -1 at 
1900 Hz (Figure 6.7), the predicted aode amplitude for the second radial aode 
(n-1) is much higher than the aaplitude for the first radial aode (n*0). The 
neasureaents , however, show the n«0 node to be slightly higher than the n-1 
node. 


In utilizing the suppression prediction prograa, it is essential to 
cheek that the total acoustic energy is conserved aeross segment interfaces. 
The nodal analysis requires the acoustic pressure and velocity at all radial 
locations to be continuous aeross segment interfaces. However, the pressures 
(and consequently the velocities) on the two sides of the interface are 
represented as expansions in sets of eigenfunctions which can be completely 
different on either side. If the set of eigenfunctions used is not 
nathenatically eonplete, then the expansions, the pressure matching, and 
consequently the acoustic energy matching, can be in error. For the first 
three conditions listed in Table 6.1 the predicted acoustic energy across the 
segment interfaces in the rotor 55 duct was conserved within 0.1 dB. For the 
a - 7 mode at 1900 Hz, however, an energy mismatch of 2.5 dB was observed at 
the downstream segment interface. This mismatch was reduced to 1.2 dB and 1.0 
dB by including 6 and 8 radial modes respectively in the analysis. Assessment 
of the effects of including radial aodes of order higher than this was 
prevented by limitations in the size of the computer prograa. 
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Figure 6.5. Measured and Predicted Coaplex Mode Amplitudes at the Upstream Plane - 

Backward Propagating Modes - 1000 Hs, m » -1, M - .21, Uniform Flow, Soft Wall 
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Figure 6.6. Measured and Predicted Coaplex Mode Amplitudes at the Upstream Plane - 

Backward Propagating Modes - 1500 Ht| a ■ -1, H ■ .32* Uniform Flow, Soft Wall 
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Figure 6.7. Measured and Predicted Conplex Mode Anplltudca at the Upstreaa Plane - 

Backward Propagating Modes - 1900 Hz, o ■ -1, M * .4, Uniform Flow, Soft Hall 
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Figure 6.8. Heaaured and Predicted Complex Mode Aaplltudea at the Upatreaa Plane - 

Backward Propagating Modea, 1900 Ha, m • 7, M ■ 0.4. Uniform Plow, Soft Hall 
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Figure 6.9 shows the predieted mode coefficients (for both forward end 
becfcward propagating modes) at the end planes of the hsrdwall and treated 
segments for m • -1 aodes at 1000 Hz. At this frequency only the lowest order 
redial mode (n*0) is eut-on in the hsrdwall duet. Consequently the forward 
propagating higher order aodes (n > 0) at the "source" plane i.e. , plane 1, 
are suppressed to eery low values at plane 2 just upstream of the treated 
segment. Downstream of the hsrdwall/treatment interface i.e., at plana 3, the 
aode coefficient distributions shown are required to aateh the energy at plane 
2. Hotiee that a combination of forward nodes n • 0, 1, 2 and backward nodes 
n ■ 0, 1 at plane 3 in the treated segment is able to natch the energy carried 
by the combination of forward aode n ■ 0 and backward nodes; n - 0 and n « 2 at 
plane 2 in the hardwall segment. Thus a redistribution of aode coefficients 
takes place at the segment interface. This is due to the fact that the radial 
profiles of acoustic pressure and velocity at the interface are expressed as a 
summation of hardwall eigenfunctions at plane 2 and as a summation of softwall 
eigenfunctions at plane 3. A similar mode coefficient redistribution takes 
place at the downstream treatment/hardwall interface. 

For m » 7 at 1900 Hz most of the incident energy is contained in only 
one radial aode (n « 0) that is close to cut-off (see Figure 6.10). This mode 
has a cut-off ratio equal to 1.1. Upon incidence to the upstream 
hsrdwall/treatment interface this mode is reflected into several higher order 
backward traveling nodes. The backward modes at plane 2 in the hardwall can 
be identified as due to reflection of forward modes at the same plane and not 
as redistribution of backward modes at plane 3 in the treated segment because 
the latter have near zero amplitudes as shown. Thus the acoustic energy 
carried by the combination of the n * 0 forward mode and the backward modes of 
several higher orders in the hardwall segment must be matched by distribution 
of energy mostly in the forward modes in the treated segment. The backward 
modes at plane 3 are small in amplitude because of the high suppression in the 
treatment. At the downstresm treatment/hardwall interface the reflection of 
the n ■ 0 forward mode in the treated segment into several higher order 
backward traveling modes can be observed. 
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Figure 6.9. Predicted Mode Distributions at Segment End Planes - 1000 Hz, 


m 
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It is fait that this combination of high reflection tad suppression of 
tha naar cut-off aoda causas tha problem of energy mismatch at aagaant 
intarfacas. The mismatch is reduced by including higher order radial modes in 
the analysis. Since these modes are "cut-off* they are highly suppressed even 
in tha hardtmll segments of the duet. Their contribution to the total 
acoustic energy at the "source" and "termination" planes is therefore 
negligible and the in-duet suppression based on the forward energies at these 
planes does not change much. 

g.3 tmssusLsaim rww 

As mentioned in Section 5.2.3, the mean flow in the Ho tor 55 exhaust 
duet was far from uniform. In addition to the radial variation of velocity 
due to the thick boundary layers at the duet walls, a periodic circumferential 
variation of the mean velocity due to the persistent stator wakes was also 
observed. It is not possible to account for the circumferential variations in 
the prediction of suppression without considerably modifying the analysis. An 
attanpt was, however, made to account for the thick boundary layers by using 
tha analysis described earlier for thin boundary layers. Sheared flow 
eigenvalues were determined by using the measured velocity profile in the 
integration of the governing differential equation across the duct annulus. 

The eigenvalues and the acoustic pressure mode shapes obtained during this 
integration process were used to obtain the sheared flow mode coefficients 
from the mode probe measurements. These mode coefficients were used in the 
suppression prediction program to obtain estimates of the suppressions due to 
the treatment in the duct. 

Table 6.2 compares the measured in-duct suppressions with the 
suppressions predicted using the sheared flow mode coefficients for the sane 
conditions as in Table 6.1. For the m ■ -1 mode at 1500 and 1900 Hz the two 
suppression values can be seen to be in good agreement. Furthermore, the 
measured values differ from the corresponding values in Table 6.1 by at most 
3.2 dB. For the spinning mode a a 7 at 1900 Hz the in-duet measured 
suppression based on the sheered flow mode coefficients is very similar to the 
measured suppression based on the uniform flow mode coefficients. The 
predicted end the measured suppressions based on the sheered flow mode 
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Table 6.2. Comparison of Measured and Predicted In-Duct Suppressions 
Based on Sheared Flow Mode Coefficients 
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coefficients for this esse, however, do not cospere as well as those based on 
the uniform flow mode coefficients. This is due to the large mismatch of 
energy in the prediction routine at segment interfaces for the sheared flow 
case. In fact even for the eases where the measurement and prediction agree 
Ci.e. , m ■ -1 mode at 1500 Hz and 1900 Hz), the energy conservation at the 
segment interfaces is off by up to 0.8 dB. It is believed that this is due to 
the assumption in the analysis that the complex acoustic pressure profile in 
the duct in the presence of non-uniform flow can be expanded in t^nss of the 
uniform flow eigenfunctions. This assumption is valid only for thin boundary 
layers and can introduce large errors in highly sheared flows, furthermore 
the expression for acoustic intensity (equation 4.39) is valid only for 
unifora flow* 

In an attempt to improve the energy conservation at segment interfaces 
in the presence of sheared flow, the program was modified so that the acoustic 
pressure mode shapes generated numerically during the determination of sheared 
flow eigenvalues are used in the modal analysis instead of the uniform flow 
eigenfunctions. This modification did not introduce a consistent improvement 
in the energy conservation at segment interfaces. It was, therefore, not 
pursued further. 

6.4 SCTSITIYITT OF SUPPRESSIOl TO TREATMENT IMPEDANCE 

All the predicted suppressions presented so far are based on the 
average normal acoustic impedance of the treatment panels measured at 140 dB 
in the absence of flow (see Section 5.3.3). Treatment impedance is known to 
change due to the presence of flow (and boundary laysr) ovar the treatment and 
also with changes in sound pressure level in the duct (see References 19 and 
4). Rasults of the experimental attempts to assass the flow and SPL effect on 
the acoustic impedance of the treatment used in Rotor-55 exhaust duct are 
described in Section 5.3.2. The average SPL in the Rotor-55 exhaust duct was 
in the range 125-135 dB. This will cause the resistance of the treatment to 
be lowered relative to the value at 140 dB. However, the grazing flow over 
the treatment will cause the resistance to increase. The treatment reactance 
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was observed to become less negative due to the presence of flow. It was 
estimated that the treatment resistance and reactanee under the test 
conditions will be within 0.1 pc of the average normal impedance values. 

Table 6.3 (for uniform flow condition assumptions) shows that by varying the 
treatment impedance within this limit the predicted in-duet suppression can be 
brought into excellent agreement with the measured suppression. 

The sensitivity of the predicted suppression to the treatment impedance 
was established by changing the resistance and reactance over a range of 
values close to the averaged awasured value. The results for m « -1 mode in 
the presence of sheared flow are plotted in Figures 6.11, 6.12 and 6.13 for f 
• 1000, 1500 and 1900 Hz respectively. At 1000 Hz the suppression is almost 
equally sensitive to changes in resistance and reactance, varying by about 1 
dB over 0.4 pc. At 1900 Hz the predicted suppression is extremely sensitive 
to the reactance and even a change of 0.1 pc in reactance can change the 
suppression by 4 dB. The sensitivity to resistance is not that high, but 
greater than the sensitivity at 1000 and 1500 Hz. In Figure 6.13 there is 
also an indication of a local peak in the suppression curve at 2 • 0.45 pc. 

6.5 FAB FIELD SPPPBgSSIOB 

The acoustic field radiated from the exhaust duct was measured in the 
far field at all three rotor speeds for the two eases: a) all hard wall 

segments in the duct and b) one treatment segment in the duct. By integrating 
the SPL directivities (at the blade passage frequency) for the treated duct 
and the hardwall duet and taking the difference we can obtain a power level 
(PVL) suppression due to the treatment. Table 6.4 lists the power level 
suppressions and the uniform-flow in-duet suppressions for all three ** 
frequencies. At 1000 Hz and 1500 Hz, m « -1 is the only rotor-stator 
interaction mode generated in the duct. At these frequencies significant 
energy is carried by the (-1,0) mode. This mode is well cut-on at both 
frequencies and will have a peak lobe of radiation in a direction close to the 
duct axis. The jet flow exhausting from the duct termination (i.e., nozzle), 
however, will introduce refraction effects and cause the peak lobe to move 
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Table 6.3. Sensitivity of Predieted Suppressions To Trestaent Impedance 

Uniform Flow 
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Plgure 6.11. Sensitivity of Suppression to Treatment Resistance and Reactance - 1000 
He; m - -1, n ■ 0,1, 2,3; M - .21 Sheared Flow 
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Figure 6.12. Sensitivity of Suppression to Treatnen 
0,1,2, 3; M - .32; Sheared Flow 
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Figure 6.13. Sensitivity of' Suppression to Trestaent Resistance and Reactance - 1900 
He; a ■ -1, n ■ 0,1,2, 3; M ■ .4; Sheared Plow 



Table 6.4. Comparison of In-Duct Suppressions With PWL 
Suppressions Prom Par-Field Data 
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mjr from the duet axis. Far-fiald measurements in tba sector covering 20* 
through 110* should be able to capture the peak lobe of radiation and PVL 
suppressions based on far-field radiation pattarn should be representative of 
the in-duet suppression due to the treataent (sawing the reflection at the 
nozzle to be of the saae order for both hardwall and treataent eases) . In 
Table 0.4, the PUL suppression (based on far-field data) can be seen to be of 
the sasa order as the in-duet suppression at 1000 and 1500 Hz. 

At 1900 Hz, the interaction nodes a • 0-1 and a • 7 coexist, with the 
latter c arryi ng a greater part of the energy. A good percent of this energy 
is, however, refleeted baek into the duet at the nozzle termination. 
Measurements in the downstream hardwall section showed the baekward energy in 
the m»7 modes to be only 4 dB below the forward energy. For the a - -1 modes 
this difference was found to be nearly 10 dB. The (7,0) node has a cut-off 
ratio near 1.1 and a peak lobe of radiation almost normal to the duct axis. 

The refractive effects of the jet will move the peak lobe to even higher 
angles. This can cause the far-field. measurement field (20* thru 110*) to 
miss part of the radiated energy. The a ■ -1 mode will, as in the case of 
lower frequencies, radiate at a smaller angle to the duet axis. The 
multi-modal PUL suppression (from far-field data) is 10.3 dB compared to the 
measured in-duet suppression values of 16.8 and 22.5 dB for the a • -1 and m - 
7 modes respectively. Further examination of the causes for the variations in 
these numbers is required. 

6.6 vnrn.T vrrLtCTLOW CHARACTERISTICS 

In Section 4.1.3 a method for measuring the reflection matrix of the 
nozzle (B^l was outlined. It was pointed out that all the elements of the 
reflection matrix (including the off-diagonal terms representing the mode 
scattering) can be determined by measuring the forward and backward mode 
distributions just upstream of the nozzle for different distributions of the 
incident modes. For the rotor-55 exhaust duct, in-duct measurements were made 
with and without the treatment segments in the duets. The mode distribution 
incident to the nozzle for the all hardwall segment configuration is different 
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froa the mode distribution for tho ooo troatod segment configuration. Mods 
coefficients measured in these two eases can be used to solve equation 4.41 
sad determine four eleaents (namely 1^, * 01 X^ 0 and X^) for tbe 
nozsle reflection natrix. .The first subscript of X represents tbe radial 
node order of the reflected a»de and tbe second subscript refers to the 
radial node order of the incident node. 

. Table 6.5 lists the nozzle reflection natrix eleaents based on the 
measured node coefficients in the rotor-55 exhaust duet. Wo dear trend is 
evident from these nuabers. For a « 7 aodes at 1900 Hz, the reflection natrix 
eleaent X^ can be seen to be very high. This points to tbe strong n » 0 
backward node near the nozzle for these conditions (Figure 6.10). 

Theoretical prediction of the suppression due to treatment in the 
exhaust duct required the inclusion of higher order radial aodes (at least 
four) in the analysis. This resulted in a nozzle reflection matrix with at 
least 16 eleaents. The larger matrix cannot be ca l culated from the available 
aeasureaents. The nozzle reflection effects were therefore taken into account 
by setting CXjl ■ 0 and specifying the backward node coefficient vector for 
the "termination" vector {Qj} . 

6.7 MODE SHAPES AXD EOPT VAT.CTT JJ fggDAgCE 

The predicted in-duet suppressions using the aeasured mode coefficients 
based on sheared flow eigenvalues and eigenfunctions were compared with tbe 
aeasured suppressions in Section 6.3. Tbe calculations on which these 
predictions were based showe d discontinuity of total acoustic energy across 
segment interfaces. It is believed that this is due to the assumption in tbe 
analysis that the coaplex acoustic pressure profile in the duct in presence of 
non-unifora flow can be expanded in terms of the uniform flow eigenfunctions. 

Figure 6.14 shows a comparison of the acoustic pressure profile in tbe 
presence of sheared flow in the hardsell segment of the duet with tbe uniform 
flow eigenfunction shape based on tbe sheared flow eigenvalue. This and tbe 


146 








f • MOO Hi 

■ ■ 7 

n - 0 (FORWARD) 
IV* 0.37 
0.02 
HAROUALL 



Uniform Flow Eigenfunction 

O Acoustic Pressure Profile 
Fron Hunerical Integration 
"Actual" 


Coaparlson of the Uniform Flow Eigenfunction Vith the Actual 
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following five figures (1 .«m Figures 6.14 through 6.19) ars for the a » 7 
aodes et 1900 Hs. The flow through the duet le assumed to have boundary 
layers with 2 percent thickness and one-seventh power law velocity profile. 
The uniform flow eigenfunction for the hardwsll duet can be seen to be in 
excellent agreement with the "actual" pressure profile (obtained by numerical 
integration). 

For the lowest order aode (n-0) propagating backward in the hardwsll 
segaent (see Figure 6.15)» the unifora flow eigenfunetion accurately 
represents the complex pressure profile in a region of the duet annulus away 
froa the walls. leer the inner wall it overprediets the actual acoustic 
pressure while at the outer wall it underprediets. The differences are, 
however, assail. 

In contrast to the hardwsll segment cases, the uniform flow 
eigenfunctions in a treated segment are considerably different froa the 
"actual** acoustic pressure profiles for both forward and backward propagating 
■odes (see Figures 6.16 and 6.17). For the forward aode there is good 
similarity between the two pressure profiles in the outer half portion of the 
duet annulus. But near the inner wall even the shapes of the two pressure 
profiles are different. 

These figures show that the uniform flow eigenfunctions can be used to 
represent the acoustic pressure profile in the hardwsll segment of a duct in 
the presence of thin boundary layers at the walls. In treated segments with 
sheared flow, however, the uniform flow eigenfunetion representation of the 
acoustic pressure profile is not accurate. 

In an attempt to improve the representation of acoustic pressure 

profile (in ducts carrying sheared flow) by uniform flow eigenfunctions the 

concept of treating the boundary layer as equivalent wall impedance was 

studied. As discussed in Section 4.1.1 the radial dependence of acoustic 

pressure (p ) in an annular duct carrying uniform flow is 
r 

P <r> - J(kr) + CY(kr) (6.1) 

p ■ r iq a r 
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Figure 6.15. Comparison of the Uniform Flow Eigenfunction With the Actual 

Acoustic Mode Shape In Hardvall Duct - f ■ 1900 Hr, a • 7 n « 
0, - .37, 6* ■ .02, Backward Propagation f 
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■ • 7 
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Figure 6.19. Coaparison of the Unifora Flow Eigenfunction Based on the 

Equivalent Impedance With the Actual Mode Shape in the Treated 
Duct - 1900 Hr, m - 7, n - 0, M - .37, i* - .02, NEX - 7, c 
■ (.51 - .551), Backward Propagation 
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and satisfies the admittance boundary conditions at the walls 


and 


a 

dr 


*P t 

ar 


iltp (1 -kM) p 


r»r. 


( 6 . 2 ) 


r*r. 


- ike <1 -kM)P. 


r»r„ 


(6.3) 


r»r„ 


These equations can be combined (by eliminating C^) and rearranged in the 
fora. 


^OkHr + A 2 (BkH) + A 3 - 0 


( 6 . 4 ) 


where H is the duet height and A^, A 2 and A 3 are functions of the 
eigenvalue, axial propagation constant, duct geometry and flow Kach number. 
The eigenvalues calculated for the boundary layer cases can be used to 
evaluate the coefficient A^, A 2 and A^ and equation (6.4) can then be 
solved to obtsJn the equivalent uniform flow admittance, i.e. , 


-A,± 

eq * 2AJVH 


(6.5) 


This value of the equivalent uniform flow admittance can then be used with the 
uniform flow eigenvalues corresponding to it (or equivalently with the sheared 
flow eigenralues corresponding to the actual wall admittance) and 


c - ig^tayo yyy - 

(m - ie^A) V k p r 2 ) - (V 2 , Vl <k r r 2 ) 


( 6 . 6 ) 


with A » (1-*H) 

in equation (6.1) to obtain the radial acoustic pressure profile. 
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Figure 6.18 and 6.19 show comparisons of the actual mode shapes with 
the unifona flow eigenfunctions based on the calculated equivalent ispedance. 
For the forward node in the treated segnent the unifora flow eigenfunction 
based on the equivalent impedance is in good agreement with the acoustic 
pressure profile. This is in contrast to the comparison of Figure 6.16. For 
the backward node, however, the use of equivalent ispedance does not improve 
the representation of the acoustic pressure profile bp the uniform flow 
eigenfunction. Further investigation is needed to seek ways of improving the 
eigenfunction representation. 


6.8 


OF AFT-OTCT 




OH FAR BOISE CTHERATIOH 


Recent analysis [Reference 20] on the generation of fan tones has 
suggested that the fan noise characteristics (both node distribution and 
acoustic power) could change when the wall impedance of any section of the 
duct (at or away froa the fan plane) is changed. In the present study the 
source characteristics were determined from neasurenents in the duct. For the 
prediction of the in-duet suppression the source characteristics obtained froa 
measurements with the treatment sections in place were used. The effects of 
the duet treatment on the fan noise characteristics are thus taken into 
account in the predicted suppressions . 


By comparing the mode coefficients at the upstream measurement oiane 
for the cases of the hardwall duct and the duct with the treatment sections, 
we should be able to assess the effect of the duct treatment on the fan noise 
generation. The measured mode distributions of the fan noise in the hardwall 
duet and in the duet with a treatment section (downstream of the fan plane) 
are shown in Figures 6.20 thrugh 6.23 for the four test conditions of Table 
6.1. The mode coefficients of the lowest order radial mode (n-O) for the 
m ■ -1 modes (Figures 6.20, 6.21 and 6.22) in the case of the hardwall duct 
differ by at most 2 dB from the mode coefficients in the presence of the 
treatment in the duct. This suggests that the treatment in the exhaust duct 
of the fan does not have a significant effect on the fan noise generation. 

For the (7,0) mode (Figure 6.23) the difference in mode coefficients is of the 
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Figure 6.20. Complex Mode Amplitudes sod Phsse Values st Upstream Plane * 
lOOO Hi, « ■ -1, M • .21, Uniform Flow, Forward Propagation 
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Figure 6.22. Coaplex Mode Aaplitudes and Phase Values at the Opstreaa Plane 
1900 Hz, a “ -1, M “ 0.4, Uniform Flow, Forward Propagation 
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order of 5 dB. The differences in node coefficients for the (7,1) cod 7,2) 
modes ere also of the sea* order thus indicating e change in the a&bient or 
fan operating conditions. 

It should be Mentioned that the data for the hard wall duct and the duct 
with the treatasnt section were obtained on different daps and the operating 
conditions were not identical (except for the fan rpa). The comparisons of 
the node coefficients in Figures 6.20 through 6.23 are thus appropriate for 
qualitative analysis only. Votiee that the treatment in the duct was located 
a few duet heights downstream of the fan OCT plane and therefore say not have 
a significant effect on the fan noise generation. 
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7.0 sssKissm mim 


Based on the theory-data comparisons presented above it can be 
concluded that the modal analysis for sound propagation in segmented annular 
duets can be used to predict the in-duct suppression due to treatment in the 
exhaust duet of a turbomachine. The modal distribution of the acoustie source 
and the reflection characteristics of the duet termination are required for 
the prediction. These ean be obtained from measurements (at least in a 
laboratory duet) or from analytical methods. 

The predicted values of the in-duet suppressions based on measured 
source characteristics and treatment acoustie impedance were found to be in 
good agreement with the measured suppressions. This is true in spite of the 
presence of thick boundary layers and circumferential flow variation in the 
test duet, two effects not fully accounted for in the theory. The presence of 
thick boundary layers did, however, cause numerical problems in the 
calculations manifested as energy mismatches at duct segment interfaces. This 
problem will require further investigation. The predicted suppression was 
found to be quite sensitive to treatment impedance indicating the importance 
of accurately determining this parameter. 

The analysis should include at least two radial nodes above those 
theoretically "cut-on" in the hardwall segment of the duct. This permits a 
valid expansion of the acoustie pressure in terms of the mode eigenfunctions 
and maintains conservation of acoustic energy across segment interfaces. More 
modes say be required when a significant amount of energy is carried by a mode 
that is close to cut-off. 


Pft£C£DB*G 


2LANX fw : 


PR£C£DL\C iV-.t;:; 


NW I'iLUyZ) 


163 


This Page Intentionally Left Blank 



wasssssi 

1. Zorumski, U. , "Acoustic Theory of Axisyametrie Hultiseetional Duets,” 
■ASA TR R-419, May 1974. 

2. Kraft, B.B. , "Theory and Beasureaent of Acoustic Wave Propagation in 
Multi-Segaented Rectangular Flow Duets," Fh.D Dissertation, University of 
Cincinnati, June 1974. 

3. Joshi, M.C., Kraft, K.B. , end Son, S.Y. , "Analysis of Sound Propagation 
in Annular Duets With Segmented Treatment and Sheared Flow," AXAA Paper 
80-0123, January 1982. 

4. Dean, D.P., "An In-Situ Method of W*. ' Impedance Beasureaent in Flow 
Duets, Journal of Sound and Vibration, Vol. 34(1), pp. 97-130, 1974. 

5. Kooi, J.W. and Sarin, S.L., "An Experimental Study of the Acoustic 
Zapedance of Helaboltz Besonator Arrays Under a Turbulent Boundary 
Layer," AIAA Paper 81-1998, October 1981. 

6. Zandbergen, T. , "On the Practical Use of Three-Microphone Technique for 
In-Situ Acoustic Impedance Measurements on Double Layer Flow Duct 
Liners", AIAA Paper 81-2000, October 1981. 

7. Syed, Asif A, "Acoustic Impedance Measurement, Patent Application, Ho. 
499, 441 filed on Hay 31, 1983. 

8. Joshi, M.C. and Kraft, K.E., "Turbofan Aft Duet Suppressor Study - 
Program Listing and User's Guide." 

9. Fiske, 6.H. , Syed, A.A. , Motsinger, R.E., Joshi, M.C. , Kraft, R.B. , 
"Contractors Data Report of Mode Probe Signal Data', MASA Contract Ho. 
■AS3-22744 . 


"CEDING PAGE BLANK NOT 


165 


10. Kraft. R.E., and Walla, W.S., "Adjointness Properties of Diffaraatial 
Syateaa with Eigenvalue Dependant Boundary Conditions With Application to 
Flow Duet Acoustics,** J. Acoustical Society of Aaarica, Vol. 61, lb. 4, 
April 1977. 

11. Cantrell, R.B. and Hart, R.H. , "Interaction Between Sound and Flow in 
Acoustic Cavities: Hass, Momentum and Bnergy Considerations," J. 
Acoustical Sociaty of Aaarica, Vol. 36, lo. 4, 1964, pp. 694-706. 

12. Zandbergen, T. , Logan, J.l. , Zaaaans, H.J., and Sarin, S.L., "In-Flight 
Acoustic Hessuraasnts in the engine Intake of a Fokker F28 Aircraft", 

ATAA Paper lo. 83-0677, April 1983. 

13. Moore, C.J., "In-Duct Investigation of Subsonic Fan 'Rotor Alone* loise," 
Journal of Acoustical Society of Aaarica . Vol. 51, lo. 5 (Part 1), 1972. 

14. Kline, S.J. and HeClintock, "The Description of Uncertainties in Single 
Saaple Experiments," Mechanical Enaineerlnx. January, 1958. 

15. Moore, C.J., "Heasureaents of Radial and Circumferential Modes in Annular 
and Circular Fan Duets," Journal of Sound and Vibration . Vol. 62(2), pp. 
235-256, 1979. 

16. Morse, M.M. and Ingard, K.U., Theoretical Acoustics . McGraw-Hill, Sew 
York, 1968. 

17. Taylor, K. , "A Transformation of the Acoustic Equation With Iaplieations 
for Wind-Tunnel and Low-Speed Flight Tests", Proceedings of the Royal 
Society, 363, 271-281, 1978. 

18. Tyler, J.H. and Sofrin, T.C., "Axial Flow Compressor Voise Studies," SAE 
Trans . . 70, pp. 309-332, 1962. 

19. Hersh, A.S. and Walker, B., "The Acoustic Behavior of Helaboltz 
Resonators Exposed to High Speed Grazing Flows," AIAA Paper 76-563, July 
1976. 


166 


151 



20. Hungur, P. and Gliebe, P.I. , "Acoustic Analysis of Fan Tons Generation 
and Propagation" , AIAA Paper 82-0769, A2AA 8th Aeroacousties Conference, 
Atlanta, GA, April 1983. 

21. Fiske, G.H. , Syed, A. A., Hotsinger, H.E., Joshi, H.C. and Kraft, B.E. , 
"Measurement of Acoustic Modes and Vail Taped snee in a Turbofan Exhaust 
Duet", AIAA Paper Mo. 83-0733, April 1983. 

22. Joshi, M.C. , Kraft, S.B. , Fiske, G.R., Sped, A.A. and Hotsinger, S.B. , 
"Sound Propagation in Seyrmted Exhaust Duets - Theoretical Prediction 
and Comparison with Measure: lentss", AIAA Paper lo. 83-0734, April 1983. 


167 ‘ 



This Page Intentionally Left Blank 



This Page Intentionally Left Blank 



AFPnOXX A 

SCATTER XI MODAL COSFF1CXEBT FILES 


Extensive work was done to determine the reliability of a particular 
probe through the elialnsion of that probe in the nodal decomposition scheae 
and its subsequent effect on the results. This could be done by eliminating 
one or two probes from the input data to the modal decomposition program and 
solving for less radial mode orders. This was done for all probes and then 
the modal coefficient files were examined for scatter. To determine the 
degree of scatter, all modal coefficient files that were generated without a 
particular probe were combined and each mode had its average value and 
standard deviation calculated. The reasoning was that if a probe was not 
functioning properly its elimination from the data would result in reduced 
scatter (lower standard deviation) assuming all the other probes ware 
functioning well. In the ease where a probe was calibrated to read high as 
was the ease in the 7600 rpm hardwall run, this proved to be an excellent node 
of detecting a bad probe. However, in the case where a probe was slightly 
niscalibrated or reading low, this method was not as effective. For example, 
in the ease where a broken diaphragm gave no reading, as was the case in all 
downstream treatment runs, the scatter was less by leaving out this probe but 
only slightly less. Figure A1 shows a comparison of scatter versus probe 
elimination for these two eases. Since the degree of confidence in using this 
method was low for low rsading probes, it was not utilized in these cases. 

In summary, modal decomposition was carried out using all functioning 
probes, except in the case of the 7600 rpm hardwall run. In that ease, the 
cause of misealibration was known before the scatter techniques were applied, 
and the scatter techniques confirmed that fact. In the downstream treatment 
eases the probe with the defective transducer was eliminated from modal 
decomposition. That is, the scatter criteria was not used in the data 
reduction process. 
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Figure Al* Scatter in the Modal Data Due to Probe Elimination at the Downstream 
Measurement Location 



APPKHDIX B 

AH ASSESSMBHT OF EBBORS ZH THE H-SITU HEASUBEHEHTS 
OF THE ACOUSTIC DffEDAHCE 


'T"T 


•T*cm Sb««t 
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P A (2) Slpil at A 

p_ (2) Signal at B 
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^Acoustic Transducers 


The sketch shore shows two transducers sens ins signals at the surface 
of the porous face sheet (A) and at the backvall of the cavity (B) . 

The Acoustic Impedance of the locally reacting SDOF panel in terns of 
the two signals 9 its thickness d and the frequency f is given by 


C(f> - -i exp {i ^(f)} cosec(kd) (Bl) 


where k - 2«f/c, c is the speed of sound and 


exp {i - P A (f ) pJ<« / P B «> P B Cf) 


The transfer function H^Cf) and the phase difference are 

obtained from the analysis of the two signals; cosec (dk) is computed from 
specified values of the panel thickness d v and temperature T in the cavity. 
First consider the error in the measured impedance due to errors in the values 
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of d and T. It can be shown that an error in the Impedance value if, (f ) 
due to an error i (ltd) is given by 

«C x (f) - 

• -(R+iX) cot(led) . «<kd) (B3) 

where 


«Ckd> - 0.000313 f <B4) 

T is in *1 and d is in millimeters. 

It is clear from (B3) that the error 4f^(f) Is zero when cot(kd) is 
zero, that is, close to the tuning frequency of the SDOF panel. Saaple 
calculations of the expected errors in the measured impedance of an SDOF panel 
due to errors in the specified thickness and te^erature are presented in 
Table Bl. The errors in the resistance and reactance values due to errors in 
the tesperaturu and thickness are quite large near the antiresonance frequency 
given by c/2d. Jl way from this frequency, these errors are quite small and may 
be ignored. In practice the cavity thickness d aay be measured with much 
greater accuracy than specified in Table Bl. The accuracy of temperature may 
be improved by accurate calibration of the thermocouples used in the 
laboratory. 

The other sources of error in the measured impedance are errors 

4H._(f ) and $♦,_(£) in the measurement of ths transfer function and 

the phase difference of the signals. These errors aay be either of the bias 

type or of randosi type as discussed in Reference 22. These types of errors 

may be minimized by increasing the number of averages in the data analysis as 

well as by reducing the bandwidth. However significant errors in the measured 

values of H._(f) and ♦ (f) may exist due to the following: 

AB AB 

1. Transducers not flush with the surfaces at A and B. 

2* Leakages fron the cavity other than through the porous sheet. 
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250.0 

500.0 

750.0 

1000.0 

1250.0 

1500.0 

1750.0 

2000.0 

2250.0 

2500.0 

2750.0 

3000.0 

3250.0 

3500.0 

3750.0 

4000.0 

4250.0 

4500.0 

4750.0 

5000.0 

5250.0 

5500.0 

5750.0 

6000.0 
6250.0 


- 0.698 

- 0.330 

- 0.031 

0.236 

0.495 

0.766 

1.073 

1.459 

2.019 

3.050 

6.273 

- 27.150 

- 3.213 

- 1.233 

- 0.425 

0.058 

0.411 

0.702 

0.967 

1.226 

1.501 

1.818 


0.006 

0.006 

0 . CJ 6 

0.005 

0.004 

0.003 

0.001 

J.001 

3.004 

- 0.008 

- 0.015 

- 0.024 

- 0.044 

- 0.109 

0 . 580 - 

0.088 

0.648 

0.031 

0.022 

0.014 

0.008 

0.002 

- 0.005 

- 0.013 

- 0.023 


.001 
01 
01 
01 
0.000 
0.000 
0.000 
- 0.000 
- 0.001 
- 0.001 
- 0.002 
- 0.003 
- 0.005 
- 0.013 
0.067 
0.010 
0.005 
0.004 
0.002 
0.002 
0.001 
0.000 
- 0.001 
- 0.001 
- 0.003 


- 0.706 

- 0.147 

- 0.033 

0.003 

0.015 


0.004 

- 0.016 

- 0.049 

- 0.104 


- 0.009 

- 0.004 

- 0.002 

- 0.001 

- 0.000 

- 0.000 

0.000 

- 0.000 

- 0.001 

- 0.003 

- 0.006 

- 0.014 

- 0.039 

- 0.197 

- 4.532 

- 0.081 

- 0.017 

- 0.004 

0.000 

0.002 

0.002 

0.000 

- 0.002 

- 0.006 

- 0.012 
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3. Inaccuracies in the calibration of transducers 


4. Inaccuracies due to insufficient aignal-to-noise ratio. 

Tinder conditions of grazing flow, the aeasursaent of Ispedtnce is found 
to be extr—sly sensitive to the position of transducer A relative to the 
surface of the porous face sheet, experience showed that for accurate 
■assure— at the fees sheet transducer aust be absolutely flush with it. 
Leakmges were eliminated by constructing • single cavity with hard walls. 

Inaccuracies in the calibrations of the transducers and insufficient 
signal-to-ooise ratio are the aost significant sources of error in the 
■assured iapedance. The calibrations of the transducers, the associated 
signal conditioning equipaaat and the A to D converter of the data analysis 
system are required because of significant differences in their response 
characteristics. The method used for the calibration of the transducers, the 
signal conditioning amplifiers and the analog-to-digltal converters consisted 
of placing then at the end of a tube flush with a hard wall termination and 
exciting a broadband acoustic signal at the other end. Under these 
conditions, both the transducers were subjected to plane wave signals at all 
frequencies up to the cut-on frequency of the (1,0) node. The transfer 
function and the phase difference aeasured under these conditions represented 
the calibration to be used vten ■assuring iap ed an c e. Clearly, the accuracy of 
this calibration woe 1 - depend on the signal- to- noise ratio during its 
— asuremeat. Figure B1 shows the calibration of the transducer system used in 
the In- Situ measure— nts in the erasing Flow Duct. The reliable part of this 
data lies roughly between 400 Hz and 2800 Hz. At frequencies outside this 
range, signal-to-noize ratio is e problem. The difficulty lies in poor 
coupling bet ween the speaker and the tube at these frequencies and a 
relatively high level of electronic-system noise. A way to eliminate this 
problem would be to use discrete frequency signals for calibration. 

r* 

From error analysis, it can be shown that the error in the 
measured impedance due to errors <H and !♦__ is given by 

AJB 
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Figure Rl. Typical Calibration of the Two Acoustic Transducers llaed in 
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Where X and X are the acoustic resistance and r e acta n ce values respectively of 
the panel be ins tested. Table B2 shows saaple calculations of the expected 
errors in the aeasured iapedanee due to given errors in the aeasured phase 
difference and the transfer function of the signals at 4 and B. A snail error 
in phase has little effect on the aeaSbred reactance over the entire frequency 
range but gives rise to significant errors in the a easur e d resistance at 
frequencies where the rea c tan c e is large. A sm a ll error in aeasured transfer 
function gives rise to a small error in both the measured resistance and 
reactance values. 

I 

! 

Over the frequency range of interest in this study, the errors in the 
aeasured iapedanee due to errors in the aeasured values of H (f) and 

AB 

t^glf) were considered to be insignificant. 
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Table B2. Expected Errora 1 b Measured I nped a n c* Doe to Error* 

in the Meaaured Tranafer Puactloa (Has) aad Phaae Difference (4AB) 

• R ■ 0.4 pc 

• “ 0.025 radian* (1.43*) error in phase difference 

• Mu/Bab - 0.05 




<R2* t JXjj! Errors due to 

«R2h. **2B : Ptror* due to <Bab/Hab 











AFPOTOIX c 
BIGUVALOE EQtJATIOHS 


The eoaplete fora of the eigenvalue equation (referred to in Seetion 
4.1.1) for sound propagation in a lined duet carrying uniform flow is 


F(y> 


- T Cl - C> 

- T ci - C> 


2 ♦ ifikKA ~ + (BkH) 2 A 2 | 

~ 0 - iflkHA 


a(l - C) + iBkHA 


- Y (1 - c) t 4 - 0 

where 8 * acoustic admittance of the duct walla 
T - Vl 
C - r 1 /r 2 


H * r 2 * r l 


T 1 ” J a <^> V T> - J a (T) T a <fY) 

T 2 - J . Vl (T) " J »»l (T) T a 


T 3 * J a*l *« (Y) - J » (Y) V* (CY) 


t 4 “ W T) W w - W T > 


A - (1 - J V) 2 


-H tl/l - (1 - M 2 ) 

I ueni 


e . 

k 


(C.l) 


(C.2) 


(C.3) 


(C.4) 


(C.5) 


1 - M 


r^ and are the inner and outer wall radii 
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and J ( ) and T ( ) ara Basaal function* of the first and 


--second kinds of ordar *. 

As Mentioned in Section 4.1.1, the eigenvalues can be determine^ by a 
Hevton-Baphson iteration of Equation (C.l). 

Alternately, Equation (C.l) can be recast as the differential equation 


ax.l 

dx D 

- | ^ + (1cH) 2 A 2 28 | 




- t(l ~ C)i*HA(T 2 - T 3 ) g 

J 


(C.6) 

where 





[ - T(1 . C 2 )] 



+ igfcfld - O 2 * (A ♦ ^ A) ♦ TAj 



+ 2 (BkH) 2 A(A + * A) 1 




| (l - C) 2 (fr 2 - J“) ♦ iskH(l - C) [tA _ *AO_=_£2 

■ J - (BkH) 2 A 2 | 

.x 3 ! 

j (1 - C) 2 (T 2 - m 2 ) - lBkH(l - C) 

tA ♦ mA(l - o] - 

C(81tH) 2 A 2 J 

.T 4 , 

| (1 - C) 2 m(Y - 1) + iBkHAd - f)(f - y)| 

(C.7) 

and Tj, 

Tj, Ij, and ? 4 are represented in Equations (C.4). 


Eigenvalues can be obtained by integrating Equation (C.6). 

Any variation 


of admittance (8) with respect to the independent variable x can be pre- 
scribed. For a linear variation 8 » xfl^ where 8^ is the wail admittance. 
Equation (C.6) and (C.7) reduced to 



B^kHA + 2(B |# kH) 2 A 2 x 


+ T(1 - C) iB„1tHA(T 2 - Tj) 


(C.8) 


x | W<1 g C ‘ 2 [ ^ - Y<1 + C 2 ) J + lB w kH(l - C> 2 £ ® (A + *» A) ♦ taJ 
♦ 2(B w ltH) 2 A(A + ® A) x 2 J 


| (1 - C) 2 (CT 2 ’ - J“> ♦ lB w kH(l - O [*- ~ j X - <B ( 

| (1 - C) 2 (T 2 - » 2 ) - iB w kH(l - f > £ yA + aA(l - Cl > J x - C(fJ 
| (1 - c ) 2 »(T - 1) + iBJcHAU - CMC - T)x | 


To find the softvall eigenvalues. Equation (C.8) oust he integrated from 
x » 0 to x - 1. The initial values for y (at x ■ 0) are the hard vail 
eigenvalues. 
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